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Abstract: The frequency dependence of the AC conductivity is studied in a holo¬ 
graphic model of a non-fermi liquid that is amenable to both analytical and numerical 
computation. In the regime that dissipation dominates the DC conductivity, the AC 
conductivity is described well in the IR by a Drude peak despite the absence of 
quasiparticles. In the regime where pair-production-like processes dominate the con¬ 
ductivity there is no Drude peak. A scaling tail is found for the AC conductivity that 
is independent of the charge density and momentum dissipation. Evidence is given 
that this scaling tail aA_c ~ appears generically in quantum critical holographic 
systems and the associated scaling exponent m is calculated in terms of the Lifshitz 
and conduction critical exponents. 
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1. Introduction 


Insulators and superconductors are central concepts in strongly correlated electron 
materials giving rise to high superconductivity, [|I|. The region in which the 
insulator-metal transition happens, 0, is characterized by novel features of the elec¬ 
tronic conductivity, 0. Moreover, the scaling of the AC conductivity, @] has long 
been one of the benchmark features of such materials. 

Holography has provided a new paradigm and a new arena for theoretical models 
that address physics at strong coupling. The setup is semiclassical (enforced by a 
large-N limit involved)^. It is a natural framework to describe quantum critical 
systems at zero and finite density. It is very convenient to describe conductivity, a 
major observable in condensed matter. A study of holographic ground states has 
indicated that they are very diverse in their properties, i, gig, 10,0,0, 
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Most holographic systems analyzed at finite density are translationally invariant. 
Exceptions also exist, using D-brane defects and magnetic vortices, but such 
systems have been much more difficult to analyze so far. The standard symmetry 
argument indicates that the real part of the AC conductivity will have a (5(a;) contri¬ 
bution as in a translationally invariant system a constant electric field generates an 
infinite current. The 6 function is related, because of unitarity and causality using 
the dispersion relations to a ^ pole in the imaginary part of the conductivity.^ This 
5-function is distinct from the one appearing in superfluid/superconducting phases. 

The interaction with momentum dissipation agents has been discussed in rather 
general terms in ||I7| , p!8| . When the interaction with momentum dissipating centers 
is IR irrelevant, a perturbative IR calculation can determine the scaling of the IR 
DC conductivity. When the dissipation is IR relevant, it can change the nature of 
the saddle point, turning the system into an insulator, 0-0. 

The formula obtained for the holographic DC conductivity is a sum of two con- 

( 1 . 1 ) 


tributions, [23 


^ _ ^pc I diss 

(^DC — <7dC + ^DC 


The first term, has been interpreted, [^, as a quantum-critical pair-creation 
contribution as it persists at zero charge density. For the RN black hole it is a 
constant proportional to the inverse of the bulk gauge coupling constant that counts 
the relative density of charge-carrying degrees of freedom to the neutral ones in the 
strongly-coupled plasma. More recently, it was verified in that the first term in 
(IM) is the electric conductivity in the absence of a heat current. The interpretation 
as a pair creation term is then natural since charged pairs are created with zero total 
momentum and therefore not contributing to a net matter flow. 


^Useful reviews for condensed matter physics applications can be found in [||, @ and section 2 

of §. 

^This issue was analyzed in holography in 
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The second contribution in (^) is due to the effects of dissipating momentum. 
When translation-breaking operators are irrelevant, the system is expected to be 
metallic and this term is expected to give the leading contribution to the DC con¬ 
ductivity. Then, a description of momentum relaxation in terms of the memory 
matrix formalism is appropriate and shows that the conductivity takes a Drude-like 
form, though no quasi-particle description is assumed |^. Moreover, in 


5D| it was 


shown that in several cases that were studied, the drag related part of the conduc- 
tivity had the form where Q is the charge density and Tj are diffusion rates 

of the various diffusion mechanisms at play. It is strongly suspected that this result 
is general, However, there are other features that are less clear and the recent 


work has cast doubt on this simple characterization of the two terms in the 


conductivity formula. In this work we will keep using the previous characterizations 
(pair-production and Drude or dissipative part) in lack of better ones. 

This general form of the DC conductivity was seen already in pure metric back¬ 
grounds in and was generalized to dilatonic backgrounds in |]n 
as the gauge field action is the DBI action, (O is replaced by 


In both cases. 


^DC = \l (ctSc)^ + (<^DcY 


( 1 , 2 ) 


giving results compatible with (|1 . 1|) in the regimes where pair creation or drag dif¬ 
fusion dominates the conductivity. In general, we expect a nonlinear formula that 
reflects the bulk action of the gauge field. In the probe DBI cases the momentum 
dissipation is due to the fact that charge degrees of freedom are subleading compared 
to uncharged ones. This means that there is a momentum conserving 5-function but 
its coefficient is hierarchically suppressed. 

In it was observed, based on (|1.2|) , that for running scalars other than the 
dilaton and in 2+1 boundary dimensions, the drag DC resistivity, when it domi¬ 
nates, is proportional to the electronic entropy. This seems to be a general property 
of strange metals where both the measured electronic entropy and resistivity are 
linear in temperature. This was extended in to more general cases using the 
massive graviton theory, and most importantly provided a kinetic explanation for 
the correlation suggesting a more general validity. 

The general properties of holographic conductivity were further corroborated 
recently by a careful study of the holographic current-current correlators and the 
associated properties of their poles in the complex plane, ^6, ^]. On the other 


hand the search for scaling regimes in the AC conductivity has been less studied. In 
pT[] it was suggested that generalized scaling geometries can have large-ca tails that 
are scaling and the scaling exponent was calculated in EMD hyperscaling violating 
solutions. This was further analyzed in pl[|. 


In it was claimed, based on the numerical solution of the conductivity equa¬ 
tions in the presence of a charged lattice that a scaling of the AC conductivity was 
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found with the phenomenologically relevant power aAC ~ 
not survive however more accurate numerical study, H- 


This suggestion did 


The above properties of holographic DC conductivity prompted the construction 
of holographic models with strange metal behavior, |^|, |H],[0, |^]. In particular, 
in [0 a DBI action was assumed to be describing the dynamics of charge which 
otherwise moves in a AdS-Schwarzschild black hole (in light-cone coordinates). There 
was a light-cone electric held Fy+ = E whose inverse played the role of a doping¬ 


like parameter that controls the phase diagram. The formalism of was used to 
compute the resistivity as well as the magneto-resistance in this system. 

The T + behavior of the resistivity in and the T -|- behavior of the 
inverse Hall angle, observed in at very low temperatures T < 30K, where a single 


scattering rate is present, were successfully described. The model was in accord with 


the distinct origin of the criticality at very low temperatures advertised in while 
the higher temperature, T > lOOA', scaling has different behaviors between the linear 
temperature resistivity and the quadratic temperature inverse Hall angle, signaling 
two scattering rates [^. In addition to the resistivity and inverse Hall angle, very 
good agreement was also found with experimental results of the Hall Coefficient, 
magnetoresistance and Kohler rule on various high-T^ cuprates, [^|-|^|. The model 


provided also a change of paradigm from the notion of a quantum critical point, as 
it is quantum critical at T —)■ 0 on the entire overdoped region. 

In this paper we compute the AC conductivity of the holographic strange metal 
in [^. In this model the DC conductivity depends on two parameters: a scaling 
variable t proportional to the temperature that also depends on doping and a scaling 
variable controlling the total charge density J. They are defined in ( 1^.17|) in terms of 
the parameters of the holographic theory. Note that both scaling variables depend 
on the “doping” related variable, E. In it was argued that x — Xo — E ~2 has the 
same dependence as “doping” when the model is compared to experimental data. 
In particular the E ^ 0 limit corresponds to large doping while the E ^ oo limit 
corresponds to optimal doping. Therefore, at fixed temperature T and charge density 
J+, both scaling variables t and J increase with x. 

There are two main regimes on the (t, J) plane, shown in the middle plot of 
figure 1^. The one in the upper left corner we call the quantum critical (QC) regime 
(also Pair-Production regime or Charge Conjugation Symmetric Regime (CCSR)). It 
is characterized by the fact that the DC conductivity in this regime is dominated by 
the pair pro duct ion/charge conjugation symmetric contribution. The regime in the 
lower-right corner is the Drude regime (DR). It is characterized by the fact that the 
DC conductivity in this regime is dominated by the dissipation (drag) contribution. 


We have defined a parameter q in (|2.2CI|) to distinguish between the two regimes, 
g > 1 denotes the DR while g ^ 1 denotes the QC/PP regime. As the drag contri¬ 
bution to the conductivity is proportional to charge density, it follows that at zero 
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charge density (J = 0) we are always in the QC/PP regime. 

We can vary the total charge density J to enter any of the two regimes. Inside 
each regime there are two regions with distinct behavior of the DC condnctivity 

• In the Drnde regime (g ^ 1), when t 1 the resistivity is linear in t (and 
conseqnently in the temperatnre). We call this regime the linear regime. When 
t 3> 1, the resistivity is qnadratic in t. We call this regime the quadratic regime. 

• In the QC/PP regime, (g 1), when t -C 1 the resistivity behaves as p ~ 

We call this regime, the regime I. When t S> 1 the resistivity behaves as p ~ t~^. 
We call this regime, the regime II. 


In the t —?■ 0 limit the theory has an effective Lifshitz exponent z = 2 while as 
t —)■ cxo it crosses over to an effective relativistic Lifshitz exponent, z = 1, 0. What 
we hnd in onr analysis is as follows: 


1. The characteristic temperatnre scale that controls the AC condnctivity is an 
effective temperatnre T^ff (with an associated scaling effective temperatnre 
teff dehned analogonsly). This is distinct from the system temperatnre T and 
is dne to the existence of a hierarchy of interactions. Similar effects have been 
observed in electronic systems and in holography, [^-|^ . 

The relation between tg// and t is plotted in hgnre In the small t regime 
teff ^ while at large enongh t, teff — t. As was observed in it is always 

teff > t. 


2. We dehne a generalized relaxation time r by the IR expansion of the AC con¬ 
dnctivity, 

a{uj) ^ aoc + i T oj + 0{uj‘^)) (1.3) 

In the presence of the Drnde peak, this is the conventional dehnition of an 
associated relaxation time. When there is no Drnde peak present, r is still 
well-dehned, althongh in that case the interpretation as a relaxation time is 
lost. 


We give an analytical formnla for r in (|4.14|) . It takes a simple form for large 
and small valnes of the scaling temperatnre variable t. In the regime I (see £g. 
2) we obtain 

T ^ Vt (1.4) 

while in the regime II (with t S> 1) it is set by the inverse of the temperatnre 


1 

q- _ 

t 


(1.5) 


These behaviors are shown graphically in hgnre for the varions regimes. 
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3. In the Drude regime {q 3> 1) where the dominant mechanism for the condnc- 
tivity is momentnm dissipation, there is a clear Drnde peak as seen for example 
in hgures ||. 

In the PP/QC regime it is also clear from the same hgures that no Drude peak 
is to be seen in the IR of the AC conductivity. 

4. At zero charge density (PP/QC regime) there is a scaling tail for the AC 
conductivity that behaves as 

/ OJ \ ^ TT 

kl ~ - g (1-6) 

For hnite charge density this tail survives not only in the Charge Conjugation 
Symmetric regime but also in part of the Drude regime, as seen in the various 
plots of hgure as well as the ones of hgure |[ The qualitative reason for this 
is that in the presence of the Drude peak, its tail falls off as ^ and this is faster 
than u~3. Therefore the scaling tail will eventually win over the Drude peak 
for u > Uq and the only condition that it is visible is that the UV structure of 
the theory kicks-in at Uuv 3> cuo- 


5. This scaling tail of the AC conductivity generalizes to more general scaling 
holographic geometries, as previously described in [^. The equation that 
determines the conductivity is given in (|6.9|) and the equivalent Schrodinger 
problem has a potential of the form V/// = Ri + where p is the IR charge 
density (that is proportional to the UV charge density). 


In particular, for a metric with Lifshitz exponent z, hyperscaling violation 
exponent 9 and conduction exponent C, with d spatial boundary dimensions, 
we hnd that in general 


a 


r\j 




Arg{a) ~ 


TT m 


(1.7) 


with 


m = 


^ + C ~ 2 

z 


- 1 , 


( 1 . 8 ) 


There are several constraints in the parameters of this formula that are detailed 
in section Moreover this formula is valid when the associated charge density 
does not support the IR geometry. 


6. There are some special cases of ( [I.8|) that deserve mentioning. For an AdS 2 IR 
geometry, the exponent can be obtained by the z ^ oo limit in (|1.8|) giving 


m = 0. 
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For hyperscaling violating semilocal geometries we mnst take 6 ^ oo, z ^ oo 
with I = —rj hxed and obtain 


m = 


d-2 

2 


r] + l 


- 1 


d-2 

d 


V 


(1.9) 


Finally, for the gange held conformal case^ we obtain m = 0 when d = 2. 

7. We hnd that for two spatial dimensions, negative valnes for the exponent m 
are correlated with the sign of the Lifshitz exponent ^ and the strength of the 
gange held interaction in the bnlk. Parametrizing the IR asymptotics of gange 
conpling fnnction Z as 

Z r'^ , r —>■ CX3 (1-10) 

in conformal coordinates, we obtain that zk > 0 for negative valnes of m to be 
possible. 


8. In the special case where the associated gange held seeds the IR scaling geome¬ 
try, K in eqnation (|1.1CI|) is hxed as a fnnction of z, 6 and the exponent m takes 
the valne 

?>z-2 + d-e 


m = 


- 1 


( 1 . 11 ) 


and is always positive. The special case where the IR geometry is AdS 2 can 
be obtained from the ^ oo limit of (|1 .1 1|) giving m = 2. For hyperscaling 
violating semilocal geometries, we mnst take, 9 ^ oo, z ^ oo with ® = —rj 
hxed. In this case we obtain 


m = \?> + rj\ —l = 2 + r] 


( 1 . 12 ) 


9. An important issne is whether the scaling of the AC condnctivity described 
above for the general scaling geometries is controlled by the dynamics of the 
charge density, or it is decided by the nentral system. What we hnd is as 
follows: the ehective Schrodinger potentials that controls the calcnlation of the 
condnctivity has two parts. One that is independent of charge density and one 
that is proportional to the sqnare of the charge density. In the generic case it is 
the hrst that controls the UV scaling of the AC conductivity described above. 
Only if the charge density is supporting the IR geometry, then the second part 
is of the same order as the hrst and it is the sum that controls the scaling of 
the AC conductivity. This special case is also the only one we found where 
the exponent m in ( [f.6p is always positive. In the other cases it can also be 
negative, but unitarity implies always that m > —1. 

^This corresponds to the bulk coupling constant function Z{(f)) of the gauge held asymptoting 
to a non-zero constant in the IR. 
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These findings suggest that this is a generic source of scaling tails in the AC 
conductivity in holographic systems. Moreover, in generic systems this scaling is 
expected to be independent of the mechanism of momentum dissipation. 

The structure of this paper is as follows: 

In the next section we review the holographic model for a strange metal that we 
will study in this paper. 

In section ^ we derive the linearized equations leading to the calculation of the 
AC conductivity and establish the presence of an effective temperature for the charge 
dynamics. 

In section |4.1| we analyzed the equations for the AC conductivity and derive an 
analytic formula for the generalized relaxation time. 

In section | we derive analytic formula for the scaling regime and numerically 
compute the AC conductivity for different values of the model parameters. 

In section we do a general analysis of AC conductivity scaling in generic quan¬ 
tum critical saddle points. 

Finally section contains our conclusions and outlook. 

There are several appendices that provide technical details. Appendix ^ con¬ 
tains the effective DBI action and the equations of motion. Appendix ^ contains the 
computation of the open string metric while appendix P contains a change of coordi¬ 
nates in the world-sheet induced black hole metric. Appendix derives in detail the 
equations for the fluctuations. In appendix ^ we present the computational details 
of the asymptotics of the AC DBI conductivity. Finally in appendix ^ we analyze 
theories with two charge densities and the associated scaling of the AC conductivity. 


2. A holographic Model For a Strange metal 


In this section we will review the holographic model introduced in |0. In the same ref¬ 
erence a pedestrian introduction to the holographic idea and its potential condensed 
matter applications was given. The charged matter is described by a Dirac-Born- 
Infeld (DBI) action^, in the probe limit on a fixed black-hole background. This 
background is a solution of the Einstein action 


/r — 


J + 


12 

¥ 


( 2 . 1 ) 


IbvrGs 

where R and £ are the determinant of the metric, Ricci scalar and the length scale 
of the model which controls the cosmological constant, respectively. This solution is 
the standard AdS-Schwarzschild solution, 0, 

The total action is given by the sum Ie + S, where S is the DBI action in (|2.5|) 
The bulk gravitational theory lives in 4-1-1 dimensions, while the dual boundary 


^Details on the DBI action and its origin can be found in 
freedom it describes has been analyzed in [pT[. 


The type of charge degrees of 






theory lives in 3+1 dimensions at the boundary of the bulk space-time. However 
we will consider a non-relativistic limit (equivalently a null reduction in gravity by 
taking the light-cone x"*" as time). In this limit the remaining symmetry is a z = 2 
Schrodinger symmetry, one spatial dimension is lost and the boundary dual quantum 
held theory lives in 2+1 dimensions. For details we refer the reader to the original 
reference, []^. 

A crucial ingredient of the approach in was to write the AdS-Schwarzschild 
solution in light-cone coordinates 


i=3 

= g++{dx^Y + g _ {dx~Y + 2g+_dx'^dx~ + ^ gii{dx^Y + guu{duf^ (2.2) 

i=2 


with the coordinates ordered as x^ = 
d++(“) = ^dn(M)(l - h{u)) 

guu{u) 4^3^,, 

where the blackening factor is 


(x+, x_, I/, z, n). The components read 
, g—{u) = h^gii{u){l - h{u)) 

’ ” ihi' 


(2.3) 


h = 1 — [u/uqY . 


(2.4) 


In these coordinates the boundary is located at m = 0 and the horizon aX u = uq. 
Notice that the radial coordinate u is dimensionless and to ensure the non relativistic 
scaling z = 2 we assign [b] = —1 (in mass units). Although the Ads-Schwarschild 
solution has the full conformal symmetry, in this coordinate system only the z=2 
Schrodinger symmetry is manifest. 

Another peculiarity of the system is that the coordinate x'^ will be considered 
as the held theory time. Such background interpolates between a z = 1 CFT at 
high T and a. z = 2 Lifshitz-like nonrelativistic system at T ~ 0. The black-hole 
temperature is given by b‘^7r‘^i^T^uo = 1. 

The background is a 5—dimensional space time in light-cone coordinates x^ = 
x±f. The spatial coordinate x plays a special role. From a held theory viewpoint this 
coordinate frame can be seen as an inhnite boost in the x direction. Therefore the x 
dependence disappears and the non-trivial dynamics happens in the two transverse 
spatial dimensions y, z, 0 

The charged matter degrees of freedom are described by a DBI action 

S = -Af J d^x^y-det (g + F), (2.5) 

where g represents the background metric and F the gauge held dual to the U(l) 
conserved current. It is possible to hnd a solution of the DBI equations of motion of 
the form 
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^ 4 ( 0 ) = {Ey + hj^{u))dx~^ + {b'^Ey + h_{u))dx + {b'^Ex + hy{u))dy, ( 2 . 6 ) 


In the solution above, E represents an electric field in the y direction. In appendix A 
of reference 0 is computed and shown the full form of the background fields for this 
solution. Using the work of Karch-O’Bannon, the nonlinear DC conductivity 
can be found, 0 


o'dc — 


Pul 




(2.7) 


where J+ is the charge density and u^{E, uq) is the critical point in the radial direction 
where both the numerator and denominator under the square root change sign, 


It will coincide, as we will show later-on, with the world-volume horizon. It is given 
as 


= ^ 175 ^ + 1 ]''^ - 1 


1/2 


( 2 . 8 ) 


we remind 


2P%EP 

where m*(0,Mo) = is satisfied. 

To describe the key steps of the procedure to obtain (P) and 
the reader first that in the holographic set up, the boundary value of A(o) corresponds 
to a background gauge field switched on in the dual field theory. We also note that 
after inserting 


in (^), the action will depend only on the radial derivative of 
h± and hy. Therefore, the system will have three u-independent quantities that can 
be identified with the one-point functions of the conserved quantities of the dual field 
theory, 


J± = 


5S 


Jy 


6S 


5h' 


(2.9) 


b± ~^.~y 

is interpreted as the charge density, J_ related to the particle number and Jy 
with the electric current in the //-direction. 

The quantity inside the square root of the DBI action ( p.5|) is not positive definite. 
If we demand reality of the action at all the points of the bulk space-time, several 
conditions must be satisfied. The first condition is the emergence of a surface besides 
the horizon satisfying the following equation® 


M - M*(U,Mo) = 0 , (2.10) 

the second condition relates the particle number with the charge density as follows 

J- = P-J^ , = (2.11) 


2Pu‘i 


®This hypersurface was called in the literature ’’singular shell” and afterwards was proven 
to be the location of the induced horizon of the world volume metric . 
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and the last one is Ohm’s law 


( 2 . 12 ) 


Jy — O'dcE , 

with a DC dehned as in eqnation ( ^7^ . 

The condnctivity receives two independent contribntions, therefore we will split 
it as follows 


'^DC ~ ^oi^DR + ^Qc) 5 
with the following redefinitions 


ao = 

2 

^DR - 23I2^2Jj3E 
2 _ 

2V^bEPufi ' 


(2.13) 


(2.14) 

(2.15) 

(2.16) 


The term cdr is proportional to the charge density. It has been interpreted, 
, as the term coming from the drag force acting on heavy charge carriers, in a 
Drnde-like paradigm (althongh there are no identihable qnasi-particles here). This 
interpretation was extended to more general dilaton-dependent DBI actions, O- 
The term age is present even in the absence of charge density and it is interpreted 
as a contribntion to the condnctivity dne to pair creation of charges, This 

interpretation is not withont pitfalls, bnt if hts the pair-prodnetion paradigm in more 
than one ways, |^, |^. We will call from now on the regime in which o'dr'^ erge the 
Drnde regime (DR). We will call the regime in which the pair-prodnetion condnctivity 
dominates , cdr o'qcj the Qnantnm Critical regime (QC). The reason for those 
names will be better motivated in section ^ 

We introdnee scaling variables t, J to describe the temperatnre and charge den¬ 
sity respectively 

vrf'feT 9 , X 

t = , J^= , , (2.17) 


2VbE 


{2AfbyV2{bEy 


The physical qnantities depend on these variables. The condnetivities, cdr and age 
read in term of t, J as 


^2 _ 
'^DR — 

^2 _ 
^QC - 


EA{t) 

7W)' 


A{t) = + \/l + i* 


(2.18) 

(2.19) 


The condnctivity formnla (|2.13|) has a complicated dependence on the parameter 
t. It can be however analyzed in different regions of the parameter space, depending 
on whether we are in a drag dominated regime (DR) or a pair-creation dominated 


11 
















Figure 1: Left: Comparison of the function A(t) with its low and high t behavior. 
Center: Location of the four regimes in the space of parameters (J, t) log-log scale. The 
black line q = 1, separates the DR respect to the QC regime. The magenta line represents 
the region with q = 10“^ and the red one q = 100. Right: Contour Plot of the DC 
conductivity as a function of the scaling charge density and the temperature. 


regime (QC). To characterize the later we will dehne the ratio among Drude and QC 
conductivities, 

q{t, J) = ^. (2.20) 

^QC 

These two regimes are illustrated in the second log-log plot of hgure |^. The DR and 
QC regimes can be split into two extra regions, namely 11 and t 1 respectively. 

In the hrst plot of the hgure we illustrate in a log-log plot the function A{t) ~ 1 
for t < 1 and A{t) ~ 2t^ for t > 1. The asymptotic behaviors imply that the 
resistivity p = 1/a will have a temperature dependence of the form 


P = 


(cToJ) H 

\/2(Jn 


g 3> 1, f C 1 

g S> 1, t 3> 1 

g 1, t -C 1 

g 1, f S> 1 


linear 
quadratic 
regime I 
regime II 


Drude regime (DR) 
Quantum critical (QC) 


( 2 . 21 ) 

We observe in the DR (large enough J) that conductivity is linear at low tempera¬ 
tures and becomes quadratic at higher temperatures. In general the DC conductivity 
for a background with dynamical exponent z is anc ~ , [^. The background in 


question here has a ^ = 2 effective Schrodinger symmetry in the IR which is restored 
to a z = 1 conformal symmetry in the UV. 

The theory has an extra parameter, the Electric held, E. This acts as an external 
parameter, that is similar to doping, pressure, electric and magnetic helds in strange 
metals. Varying E, from (|2.17|) we observe that this changes t and J and therefore 
ahects the crossover of the conductivity from linear to quadratic. Indeed, in |0] it 
was shown that varying oo < E < Q produces a phase diagram for the drag-related 
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conductivity that is qualitatively similar to that measured recently in cuprates, [5C 


Here E ^ oo corresponds to optimal doping while E ^ 0 corresponds to the over¬ 
doped limit. The cross-over scale between linear and quadratic regimes asymptotes 
to zero as £' —)■ 0 while it diverges when E ^ oo. 

In []^ the general magnetoresistivity was also computed. It was shown that it 
reproduces the T + T‘^ behavior of the inverse Hall angle observed in at very low 
temperatures T < 30i^, where a single scattering rate is present. 

In addition to the resistivity and inverse Hall angle, very good agreement was 
also found with experimental results of the Hall Coefficient, magnetoresistance and 
Kohler rule on various high % cuprates ^ ^ ^ 

I MM p| . This model provides a change of paradigm from the notion of a 


59 


quantum critical point, as it is quantum critical as T —)■ 0 on the entire overdoped 
region. In this sense this work breaks apart from other holographic approaches 
[p3| , |6^ |6^ , where the measured transport is due to loop fermion effects. As such, 
it is applicable to a more general class of materials e.g., d and /-electron systems, 
where the low temperature resistivity varies as T + T^ [^| and exhibit a quantum 
critical line 


3. Stationary transport from Kubo formulas and the effective 
temperature 

We will use linear response theory to study the electric properties of the system. To 
do so we need to switch-on fluctuations for the gauge field on top of the background 
discussed before A = A(o) -|- a with a inhnitesimal. These fluctuations will evolve 
with a dynamics determined by the effective Lagrangian (see appendix ^ for a 
derivation) 

Leff = —AT V— det s (, (3.1) 

8v-dets / 

where / = da, the effective (open string) metric is defined as smn = Qmn — 
{E(^o)g~^E(^o))MN, the effective coupling is g^ = \/— det s /det {g + F(o)) and the 
vector present in the Chern Simons is 

p _ det {g + K(o)) qMNqPQ /q o'! 

Vij —-3- ^MNPQR y y , [o-^) 

O 

^ [[{g + F^0))-T^ -{M^ AT)] . (3.3) 

Notice that we are raising indices with the open string metric smn- In particular 
the line element is 

+ 2s+_dx~^dx~ + 2s+udx~^du + 2s+idx~^dx^ + s _ {dx~Y -|- 

2s_idx~dx'' + Sijdx^'dx^ + Suudu^ + 2suidudx^ , (3.4) 
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where the exact form of the components of the world-volume metric are presented in 
appendix P in terms of the background metric and gauge held. Finally, the equations 
of motion for the huctnation are 

= 0, (3,5) 

3.1 The effective temperature 


We will now show that the open string metric has a horizon in a modihed location 
with respect to the background horizon located at uq. In particular we will show 
that M* represents the location of the horizon for the effective metric s. 

As a consequence of this fact, the charged matter is at a temperature different 
than the bulk system. This can happen in other holographic systems, [^, 

E 


This can happen in other holographic systems, 
and is a peculiarity of the large-iV limit implicit in holographic theories, that 
generates a hierarchy of interactions. 

To show the previous statement we will diagonalize the metric (|3.4|) . In doing 
so we will be able to extract the location of the horizon and also we will be able 
to compute the temperature of the induced black-hole. We change coordinates as 
follows 


dx"*" = dr -I- f^(u)<lX + /+ (M)dM , 
dx~ = dX" -h f~{u)dT + f~ {u)du , 
dy = dY + f^{u)dT + f^{u)du , 

d^; = dZ -|- f^{u)dT , (3.6) 


and also rescale the radial coordinate as m —)■ The full expression for the 

functions /’s can be found in appendix |^. The diagonal metric becomes 

ds^ = Suudu^ + s _(dX“)^ -1- Srrdr -|- SyydY'^ -|- SzzdZ"^ . (3.7) 


where again the component functions can be found in appendix |y. From the new 
form of the metric we can hnd the location of the open string horizon and the effective 
temperature associated to this black hole. We hnd that Srr(l) = s“,^(l) = 0, which 
indicates that the new horizon is located at the position m = 1 or equivalently at m*. 
Expanding around this point we obtain the near-horizon expressions 


(2 


Un 




4(1 - u) {2{J+yulu^ 6u^ 

2(1 — u)uq {2ul 


5ulul + ut 




bH^uU{J+yu, + uo 


(3.8) 


(3.9) 


s _~ 




Yb'^Ui, 

7^’ 

1 e 


(3.10) 

(3.11) 
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t t 

Figure 2: Dependence of the world-volume temperature T^ff with the parameters J and 
t, the effective temperature is always bigger than the background temperature T. 




where we have introduced the redefinition J_|_ = J^. Following we compute 

the effective temperature T^ff from these formulae 


^e// - ^ ^ 


' 6mo + 2{J+yulu^, — SmqM* + M* 
2ul{uo + (J+)2m*) 


(3.12) 


The previous equation written in term of the scaling variables (|2.17|) reads 


teff 



J'^A(tYA{t) — 2^2 + (lAA(t) + 2t^ + 3) 


(3.13) 


We illustrate in Fig. ^the t dependence of 4//- 

Note that the function teff{J,t) has a minimum around t ~ 1. It is possible to 
find its exact location in the limits of J —)■ 0 and J —)■ cxo 


teff{tmin) = 1 , = Oo) = 2 ^ 0.59 (3.14) 

\/\/l33 cos(0) — 1 

teffitmin) ~ 1.274 , t™„(J = 0) = 2-^^^^ 




0.85 


where 0 = | ^tt — tan“^ 

It is worth mentioning that te// > t always, as first observed for holographic 
systems in [^. Therefore the presence of the electric field heats up the system 
above the temperature of the heat bath given by the bulk black hole. 

If we analyze 4// in the four distinct regimes of the parameters J and t (defined 
in ( p.21|) ) we obtain: 


® There is also a second minimum that satisfies < tmin and —)> 0 when g —0. 
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^eff 


2-3/4^-i/ 2 linear regime 

t quadratic regime 

regime I 

t regime II 


(3.15) 


4. The fluctuation equations 


In the previous section we derived the equations of motion for fluctuations on the 
black hole background with Schrodinger symmetry in equation ( |3.5|) . We will work in 
the new coordinates introduced in (|3.6|) . As the background electric held is pointing 
in the ^/-direction we will switch-on a huctuation in the 2 ;-direction which will be 
transverse. We will solve this equation and from the source and vev of the solution 
we will determine the two point function of the dual current and from it the AC 
conductivity using the Kubo formula. 

To compute the frequency dependence of the conductivity it is only necessary to 
study the time dependence in the helds. Without loss of generality we will switch on 
a huctuation in the 2 ; direction = az{u,x~,Ty. Since in the dual held theory the 
time coordinate is given by x^, we need to take into account this fact. To do so, we 
write the near boundary expansion of the change of coordinates in Eqs. ( 3.6|) as 


(a ^ = r - -^X 

y = Y , z = z, 


X = X + p_T , 


(4.1) 


We therefore introduce the following plane wave ansatz to evolve the helds with the 
time x~^ 

f{u,X-,r)^fiu)e-^^^tirX-)_ ( 4 . 2 ) 

The only non-trivial Maxwell equation is 


—h ) a^iu) = 0 , (4.3) 


and the Chern-Simons term vanishes for this sector. In order to use the holographic 
prescription to compute the retarded Green function, we need to impose as a bound¬ 
ary conditions a^(0) = oq and to select the infalling mode at the horizon. 
Performing a near-horizon expansion of Eq. (|4.3|) we obtain 


/{AirTeffY . . 1 ^ 

— 7 -— az{u) H + ^^(w) = 0 

[u — I)"' u — 1 


(4.4) 


^In appendix ^ we write the full system of equation of motion 
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The near-horizon solutions are az{u) ~ (1 — with the dimensionless frequency 

w dehned as 


• (4.5) 

dvrTe// 

Choosing the negative sign ensures the selection of the infalling boundary condition. 
Notice that the scaling frequency is naturally dehned with T^fj. 


4.1 The DC conductivity and the relaxation time 

To solve the differential equation for arbitrary frequency it is necessary to use nu¬ 
merical tools. The low-frequency conductivity however can be obtained using per¬ 
turbation theory. To do so we will expand the held az up to second order in w and 
the infalling boundary condition will be set by expanding the appropriate solution 
of eq. ([4.4|) as follows 

az{u) = (1 - u)-^^ + wA^Piu) + w^Af\u)) . (4.6) 


Using the previous ansatz and the Kubo formula we can write the conductivity as 




2bJ\f7TT i a'^(O) 

u az{0) 

[ 47rTX^°^(0) t ( X^'^(O) V To; A^^^(0) \ 

y u ao Teff y ao J ao ) 


//n\ 

In order for the system to have a hnite DC conductivity Az "^(O) must vanish. 
We do know already that the present model has a hnite conductivity, from its cal¬ 
culation in 1^, therefore assuming .4.i°^(0) = 0 and setting ca = 0 we obtain the DC 
conductivity in terms of the solution in ([4.7|)^ 


a 


DC 

ZZ 


T { a^^\q) 

Teff Y ao 


(4.8) 


There is another transport coefficient with units of time that can be read oh eq. 
(t4.7|). We write the AC conductivity as 


cr(a;) ^ cdc (l + + C>(a;^)) , (4.9) 

1 i T 

ao 

If the conductivity is gives as a sum over poles in the upper half plane, then r 
is determined by the pole with the smallest imaginary part. In a regime which is 
of the Drude type, r can be interpreted as the relaxation time. We will call it 

®In any case we will solve the system and we will confirm that = 0. 
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the generalized relaxation time. In other regimes bit cannot be interpreted strictly 
speaking as a relaxation time. 

To hnally compute the conductivity and the generalized relaxation time it is 
necessary to calculate the boundary value of the derivative of the hrst and second 
order helds. 

Substituting the expansion ( |4.6|) into eq. ( ^4.3|) the system can be rewritten as 
follows 

du {aQ{u)a{u)duAf iu)) = , (4.11) 

where the index i takes values i = 0,1,2. is a source term depending on the solu¬ 
tions at lower order in the perturbative expansion, and with the hrst term vanishing 
Sz~^^ = 0. The a functions are dehned as 

aoiu) = ;^(-l + ) «(“) = 


The function a{u) was normalized so that a{u) | j^^.o,£;^o = 1- The source terms and 
the explicit calculation to solve the differential equation can be found in appendix 
1^. The general solution for eq. ( [4.11|) can be written as 




u 


= Cf + 


du 


ao{u)a{u) 




(4.13) 


the integration constants are hxed demanding regularity for the solution at the 
horizon and is directly identihed as the non-normalizable mode. 

Finding the solution for [4.11| and substituting into Eqs. (|4.8|) and (|4.10|) we 
obtain the expected value for the static conductivity Eq. (p.l3|) and the following 
expression for the generalized relaxation time 


Teff T = 


log 2 


iu 


3/2 


T, 


2n 


— O'DC 


eff 


diihN'Ui, T 


"dx^log 


'0 


a{xy 


1 — X 
1 -I- X 


+ 


1 Tf 


/ dxp'(x)log(l-x^), 

(OncA'nlui' T Jo 


(4.14) 


where 


p{u) = -2g^ V-s 


{-l + u‘^)uo. 
b‘^ij\fu 


P^- 


+ F 


The formula (|4.14|) gives an analytic expression for the generalized relaxation time. 
In the limit where both the charge density and the electric held asymptote to zero, 
the integrals vanish because p{u) = a{u) = 1. In this limit the generalized relaxation 
time takes the constant (QC) value 


log 2 
2ttT ■ 


(4.15) 
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r 

J = 0 

^ r-‘ 

t > 1 

t < 1 


Table 1: High and low-temperature behavior of the generalized relaxation time 


When the charge density J = 0 there are also two interesting limits we can study. 
The hrst one is t <C 1 where 


a{u) ~ 


p{u) ~ 


^ A. 

ju'^ + U^ + 2u^ -I- M -f 1 

(4.16) 


u + 1 

V2t^ 

1 u+1 

(4.17) 

1 + 2u‘^ + u + 


and substituting these functions in ( [4.14|) we can write the generalized relaxation 
time as 


Tt 


t3/2 rl ^ 

^(1 + m ) ( m 2 + 1 ) (u^ + u + l) 


^ 0.1565 


t < 1. 
(4.18) 


The other regime we can study is t S> 1 where the functions a and p read 


a{u) 

p{u) 


u{u{u + 2 ) + 2 ) + 2 
^ St\u + 1) 

^ -u(m +1)^ + 1 
8t^{u + l) 


In this case the generalized relaxation time is given by 


Tt 


log 2 

2ti 


log(32) - 2 
327r 


t-4 + (9(^-8) 


t > 1, 


(4.19) 

(4.20) 


(4.21) 


Finally we summarize in the Table |4.1| the two asymptotic values obtained from the 
general Eq. ([4.14| ) and obtained above. 


5. The AC conductivity 

To understand the full frequency dependence of the optical conductivity it is neces¬ 
sary to resort to the numerical solution of the equation. In order to have analytical 
control we can study the high-frequency behavior using WKB techniques after writ¬ 
ing the equation as a Schrodinger system. We will start this section by analytically 
studying the high-frequency regime of the system and then we will study the model 
for arbitrary frequency using numerical methods. 
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5.1 High Frequency behavior 

In order to have an idea of the behavior of the condnctivities at high frequency we 
can rewrite the equation of motion as the Schrodinger equation, from which we can 
read the effective potential that can provide insight on the solutions. To do so we 
rewrite Eq. (^4.3|) as 

Ho 


{Hia'^)' - w'^a^ = 0 , 
til 


(5.1) 


with 


Hoiu) = 


4(1 — Ci{u) 


{AnTeffY b'^i^uou p{u) 

If we redehne the radial coordinate and the held as 

dr = Ho{u)~^^‘^du 

with 

= 

the equation can be rewritten as 


Hiiu) = ao{u)a{u) , (5.2) 


az{u) = 


Hoju) 

HAuV 


—il)"{r) + {y{r) — il){r) = 0 
and the Schrodinger potential reads 


V(n) = 


Hq{u) (H[(u)4i'(u) 


</>W) 


HAu) 




(5.3) 

(5.4) 

(5.5) 

(5.6) 


where the derivatives are taken with respect to the coordinate u. 

It is not possible to do an analytic integration of Eq. ( ^.31) , however we can 
obtain the near horizon and near boundary behavior as 


r ~ — log(l — u) 

AM‘^7lTeff 0/2 

dun 


M < 1 . 


(5.7) 

(5.8) 


In the new coordinates the boundary and horizon are at 0 and oo respectively. We 
can write the asymptotic potential 


V ~ —c(t, J)e 
5 


V 


36r2 


, r > 1, 
r -C 1, 


(5.9) 

(5.10) 


where the constant c(t, J) has a complicated dependence with t and J. Nevertheless 
the formula simplihes in the QC and DR 


^ |(l + g(t, J)) 


c(t, J) ~ 


1 - 


2q{t, J) 


q{t, J) < 1 
q{t, J) > 1. 


(5.11) 

(5.12) 
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t = 0.03 , a\rag(t,J) = q a%air(t) 


t = 0.1 , a^drag(t,J) = q a%air(t) 



Figure 3: Schrodinger potential for the different regimes of interest. 


In figure ^ we show the Schrodinger potential for all the regimes of interest. 

In the QC regime, g 1, the qualitative behavior is similar in the two sectors I 
to II that correspond to f 1 and f 1 respectively (see hgure |I]). Indeed, close to 
the boundary the potential diverges as and then approaches the horizon mono- 
tonically. The behavior changes in the DR regime q ^ 1 because an intermediate 
well is formed for high enough values of q. It is worth mentioning that for g —)■ cxo 
the region in the potential is squeezed towards the boundary. Minimizing Eq. 
( b.6| ) we can compute the location and depth of the well in the limit g —)■ cxo 

V{u) ~ -u {v? — l) , g —)■ oo 

1 

Vmin — “2 X 33/2 ~ -(0.310) , 


(5.13) 

(5.14) 

(5.15) 
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V(r) 



Figure 4: In this plot the asymptotic forms of the potential are indicated. The near 
boundary region is to the left while the near horizon region is to the right. The system has 
analytic solutions in both asymptotic regions. 


with r ^ 2 (tan ^ {\/u) + tanh ^ (\/u)). From figure we observe that the location 
of the minimum satishes 

"^miniSl) ^ '^miniS^') (5.16) 


and that the depth of the well does not depend on the value of q. 

Having the asymptotic expansion of the potentials we proceed to solve the 
Schrodinger Eq. (|5.5|) analytically in these two regions (see hgure ^ and we will 
then match them assuming that the frequency is much bigger than y/1 Vmin \ ■ 

The asymptotic solutions read 


^(r) = c(f, J) (2c(f, r(l - 2iw;) , r > n (5.17) 

ijj{r) = air^^^Ji/ 3 {rw) + a 2 r^^‘^Yi/ 3 {rw) , r < tq (5.18) 


where we have hxed the integration constants in the near horizon solution selecting 
the infalling boundary condition and normalizing the solution as 'ip{r —>■ cxo) = 

The constants ai and 02 will be hxed by the matching conditions of the wave function. 
Before doing so, we will use Eq. (|5.17|) and the Kubo formula to write a formula for 
the optical conductivity in terms of the integration constants Oi and 02 


(6Wt) ^a{oj) 


/27ry/^ (y/Sai + as) /cu\-V3 

VT; r (1/3) r (7/3) 02 \ fJ 


(5.19) 


The previous asymptotic solutions are valid for arbitrary frequencies. For fre¬ 
quencies that are much bigger than the lower value of the potential inside the region 
r G (ro,ri) (see hg. ^ we can go further. In this case, the intermediate solution will 
be a combination of plane waves 


^/,(r) = die"'" + dse-^'T 


(5.20) 


22 










Matching this solution with the wave function ( 5.17|) we fix the constants 




hi = 1 


c(f, J) 


02 = 


2 w 
J) 

2 w{ 2 w + i) 


4:w{2w + i) 

2 {w + i) / ' 


(5.21) 

(5.22) 


Following the same matching process we can fix the constants oi and 02 and we 
finally obtain the analytic expression for the wave function in the near-boundary 
region. 


Oi = 


02 — 


vre “ / 


12y^ 


^Yi (x) (aie^“(l — Qix) + 6 ia 2 X + 02 ) -|- 6 xY_ 2 {x) (02 -|- 


(5.23) 


vre ( 


12y^ 


yJ^{x) (aie^“(l — 6 ix) + 6 ia 2 X + 02 ) -I- 6 xJ_ 2 [x) (02 -I- j , 


(5.24) 

with X = ro w. Our analysis is valid for although x can be either a: S> 1 

or a; <C 1 because tq can become small if q is high enough. 


We expand Eqs. (|5.23|) and (b.24|) for a; 1 hrst 


Oi 

02 


-5^x^/3r(-|)+4(fe-i)r(|) . 

24 ^ 

Svr 

~12^V^r(|) 


and obtain the AC conductivity of the form 


rota < 1, > |14im|, (5.25) 


SAfr (1/3) 
15£5/3r (7/3) 





bJklY' 

rouo J 


(5.26) 


This behavior describes an intermediate regime for very small tq. In the opposite 
case, X 3> 1 the constants asymptote to 


ai = 
0,2 = 


(- 1 ) 5/12 


(- 1 ) 11/12 



rota > 1 > \Vmin\ , 


which produce the following expression for the UV optical conductivity 


{b^iAfT) ^ a{u) ^ 


/27ry/^ (\/3 + i) /a;\-V3 

[Y) r(1/3)r(7/3) vtJ 


(5.27) 


(5.28) 


This last formula is the correct estimate for the conductivity for large enough u. 
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5.2 The numerical computation of the AC conductivity 


To compute the full frequency dependence of the conductivity it is necessary to 
find numerical solutions. In order to implement properly the infalling boundary 
condition at the horizon we must find first an approximate solution near the horizon 
a^{u) ^ al{u). 

We redefine the gauge field as a®(n) = (1 — Az{u), and using the Eq. (S) 
we find an analytic and regular solution for Az{u) in the region of u ~ 1 as an 
expansion in (u — 1) up to fourth order. 

We then integrate numerically Eq. (|4.3|) from the point u = 1 —e to the boundary 
located at the cut-off e{J,t), using as a boundary conditions a^(l — e) = 0^(1 — e) 
and 0^(1 — e) = af (1 — e). After doing the integration of the differential equation 
and using the holographic dictionary we extract the conductivities in the regimes of 
interest. In order to plot dimensionless quantities we normalized the conductivity as 
described below. 

The form of a Drude peak in the low frequency AC (cj <C tJuv) conductivity is 
given by 


uM = (5.29) 

1 — ITU 

where aoc is the DC value of the conductivity and r is the relaxation time. If 
the validity of ( b.29|) extends to a region in which tuq 3> 1 with uq uuv, the 
conductivity will obey the following power law in this regime 


a{u) 


O'DC 


U^ 


+ 


ir Vpc 
u 


—1^27r/2 


(5.30) 


In what follows we will state that the conductivity has a Drude peak if it can be 
fitted with the formula ( ^.29|) up to some frequency satisfying tuq > 1 (wq ^ (^uv)- 
On the other hand, from the UV behavior of the system under consideration, 
the conductivity (|5.28|) satisfies a power law with exponent —1/3 and phase 30°.We 


reproduce the asymptotic formula here. 


a 


27i\ 


ijz 


2{^j r(i/3)-'r(7/3)-^(^ 


^-1 


-1/3 


A/6 ^ 2 X 0.67261 ( - 
\T 


u 


n/3 




(6.31) 


with d(a;) = {b‘^£J\fT) ^ a{u). 

To infer the existence of a Drude peak in the IR, we estimate the value of 
Uo > 1 /t and we will compare it with the UV cut-off ro® (as shown in figure]^). This 
analysis is done for three different values of t in the DR and in the QC regime. To do 
so, it is necessary to compute the generalized relaxation time. In figure we plot r(t) 
for the the values J = 10"^ (DR) and J = 0 (QC). We use the analytic formula ( |4.14|) 


V 


®From figure ^ we estimate the location tq where the potential deviates from its UV behavior 
-5/(36r2). 


24 












J = 10‘' 


J=0 




Figure 5: The generalized relaxation time r as a function of the scaling temperature 
variable t.The left plot shows the behavior of the dimensionless quantity Tr (T is the 
temperature) deep in the Drude regime (DR). The right plot shows the behavior in the 
Quantum Critical regime (QCR) at zero density. Dots show numerical data and the con¬ 
tinuous line is the analytic formula obtained in the text using perturbation theory . 


and the numerical data. From the plots we observe a perfect agreement among the 
analytical data and the numerical one. We take this result as a non trivial check of 
our code. 

We now dehne ujuv ~ l/?"o and we plot it in table || together with uq for t = 
1/30,1/10, 2 in the QC and in the DR regimes^®. In the QC regime we observe that 
Uq S> ujuv, in consequence we should not expect to see a Drude peak in such regime. 
On the contrary in the DR we observe that uo <C uJuv, and a Drude peak forms in 
the present system as we shall see. 



J = 0 (QC) 

J = 10^ (DR) 

t = 1/30 
t = 1/10 
t = 2 

Wq > 1000, Wuv ~ 2 
Wq > 200, Wjjv ~ 2 
tco > 9, Wuv ~ 2 

Wq > 0.6, W^y ~ 0 
Wq > 0.3, W^y ~ 0 

Wo > 0.2, w^y ~ 0 


Table 2: wq for different values of t in the QC and DR regimes 

With the previous analysis we proceed to study the conductivities. In hgure ^ 
we plot the real (left) and imaginary (right) part of the conductivity for q = 10^ 
(DR) and for g = 0 (QC). We also htted the Drude function ( b.29|) on top of the 
numerical data. As expected from the previous analysis, in the QC regime we do 
not see the formation of a Drude peak. The point in which the conductivity deviates 
from the the Drude htting is much less than the Uq values estimated in the table 
||. We observe such deviations because the UV physics start dominating at those 

i°Note that g(l/30, lOQ - g(l/10, lOQ - and q{2, lOQ - 10^ 
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Figure 6: Conductivity for three different temperatures in the DR (up) and in the QC 
(down) regimes. The left figures show the real part of the conductivity while the right 
figures the imaginary part of the conductivity. Dots represent the numerical data, continues 
lines correspond to a fit to the Drude peak formula. 


scales. On the other hand for q = 10^ (DR) the numerical data almost agrees with 
the Drude htting in all the plotted region. That conhrms the fact that the UV scale 
happens to satisfy uuv S> cuo, and the Drude paradigm seems to be valid in this 
dragging regime even though there are no quasi-particles. 

In the following we will analyze the UV properties of the conductivities using our 
numerical results. In hgure we show in log-log plots the real (left) and imaginary 
(right) parts of the conductivity. The parameters in these plots are properly adjusted 
to display the behavior in each of four regimes of ( |2.21|) . 

• We set hrst t = 1/30 (upper plots of hgure ^ which corresponds to either 
the regime (I) or the linear regime. The values of q used for (I) are g = 0 
and q = 10“^ and for the linear regime are q = 100 and q = 1.2 x 10^. The 
conductivity was also computed for the transition value g = 1. 
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t = 0.03 , a^*,o(l,J) = q t = 0.03 , = q CT^pairW 




Figure 7: Real (left) and imaginary (right) part of the conductivity. In the top figures the 
conductivity is plotted for different values of the parameter q ranging from the QC regime 
[q <C 1) to the DR regime {q ^ 1) for a fixed value of the scaling temperature variable t 
with f <C 1. In the bottom plots q is again varied in the full range at a fixed value of t, with 
f > 1. The colored continuous lines show the Drude fitting using the analytic computation 
of the generalized relaxation time. Straight black and dashed gray lines show the UV and 
the intermediate power law behavior of the AC conductivity. 


• In the lower plots of figure 0 we show the conductivities for t = 2 and q = 0, 
q = 10“^ (II regime), g = 1 (transition) and q = 100, g = 1.2 x 10^ (quadratic 
regime). 


From hgure we observe a similar qualitative behavior in the regimes (I) and 
(II), the same happens among the linear and quadratic regime. The location of the 
UV physics in all the regimes is at ojuv ~ lO^e//- We observe in the QC the absence 
of the Drude peak and we also verify that the conductivity obeys the power law 
obtained in subsection ^T|, Eq. (|5.28|) . 

In the Drude regime we conhrm the appearance of the Drude peak in the IR. 
However, unlike the QC, we observe a transition from the peak to a power law of 
the form indicating the existence of an intermediate regime which we were 
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Figure 8: Absolute value (left) and argument (right) part of the conductivity. In the top 
figures the conductivity is plotted for different values of the parameter q ranging from the 
QC regime {q <C 1) to the DR regime {q S> 1) for a fixed value of the scaling temperature 
variable t with t <C 1. In the bottom plots q is again varied in the full range at a hxed 
value of t, with t > 1. The blue continuous line shows the Drude fit using the analytic 
computation of the generalized relaxation time. Straight black and dashed gray lines show 
the UV and the intermediate power law behavior of the AC conductivity. 


not able to extract analytically from the analysis of the Schrodinger potential. We 
note that we observe the appearance of the intermediate regime when a well in the 
Schrodinger potential is formed (see hgure so it seems to be related to it. The 
exponent —I in the DR should show up at much higher frequency, but we were not 
able to numerically verify this because at high frequencies our numerical code was 
no longer reliable. 

In Fig. H we plot the absolute value and the phase of the conductivity. The 
plots helps verify independently the presence of the power law for the intermediate 
frequency regime of a similar nature as was seen in the cuprates in [^. 

In the QC regime, g <C 1, we observe a power-law decay with an exponent 
of —1/3 and a constant of proportionality independent of t and J consistent with 
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Eq. ( 5.31 ). The phase approaches asymptotically the value of 30°. This is again 
reminiscent of the power law seen in the cuprates, 

In the DR we observe that the bigger the value of q, the closer to 90° the phase 
gets. We interpret this behavior as the consequence of the competition among the 
drag and the critical pair-production physics. As the drag contribution is dominating 
with the increasing q, the Drude description of the AC conductivity is becoming 
better. 

For q ~ 10^ (large), the conductivity behaves as u ~2 for oj > lOOTe// (linear 
regime and quadratic regime). Note that the UV behavior of the conductivity for 
such a large value for g in the linear regime is 


a~22(l + i) 


u 




eff 


and in the quadratic regime it is 


a 


8(1+ i) 


u 


Z 


eff 


(5.32) 


(5.33) 


(see Fig. ^). In view of this fact the conductivity is expected to approach a constant 
phase of 45° and this is what we actually observe in hgure This scaling happens in 
an intermediate regime and we expect that for much larger values of the frequency 
it will asymptote to the behavior that we found by the matching conditions 
method. This was verihed numerically for q = 100 and we observed the tendency 
of the data to approach the value 30° after almost a constant phase of 45° for some 
range of frequencies. 

In hgure ^ we show a check the result of the asymptotic scaling of the AC 


conductivity indicated in the analytic formula in Eq. (|5.31|) for different values of t 
and q. We have obtained a consistent result between analytic and numerics, a fact 
that provides a credibility check for our numerical calculations. 

Our analytic computation of the high-frequency behavior of the conductivities 
shows an intermediate regime in which the real part of the conductivity behaves as 
a constant when tqw -C 1 but w 3> 1 (Eq. ( |5.26| )). To reproduce such a regime it is 
necessary to go deep in the DR (see Fig. H), so we hxed J = 10^ and computed for 
different values of t both in the linear and quadratic regime. We show in Fig. the 
conductivity for t = 0.03,0.1,2,6. The upper plots are the real and imaginary part 
of the conductivity and the plots below are its absolute value and its argument. 

We observe that in the linear regime (t < 1) a plateau in the real part of the 
conductivity is formed at a; > lOOTe// and the imaginary part has a behavior as 
we showed in section |5.1| ^^. For uj < lOOZ^ff the conductivity is very well-htted by 


the Drude formula. Note for example that the phase reaches 90° as we would expect 


^^The continuous line in figure represents the Drude fit. We do know however from Eq. ( 5.30| ) 
that at high frequencies the imaginary part of the Drude form asymptotes to u!~^. 
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Figure 9: Absolute value of the conductivity in the QC (g = 10“^) and in the DR (q = 10) 
regimes for the scaling temperatures t = 0.03, t = 0.1, t = 2 and t = 10. We observe the 
same power-law, independent of the temperature and regime, in the variable uj/Tf,ff for 
high enough values of the later ratio. 


for a Drude behavior, and then presnmably decreases to 45° and stays there for some 
range of the freqnencies because of the behavior. 

The quadratic regime {t > 1) just shows the intermediate regime with conduc¬ 
tivity scaling as and phase 45°. 

In summary, we associate the presence of the Drude peak to a regime in which 
the drag related dissipation associated with strongly coupled physics dominates over 
the quantum critical pair-production mechanism for conductivity. In the QC regime 
the scaling (UV) physics (cr ~ appears at relatively low frequencies compare 

with the associated UV scales of the DR (cu > lO'^Te//). Consequently the scaling 
tail of the conductivity is controlled by the quantum critical physics. 


6. AC conductivity scaling and quantum critical saddle-points 

We will now consider the AC conductivity in holographic systems with a critical IR 
geometry characterized by a hyperscaling violating exponent 6 and Lifshitz exponent 


Till , in the absence of momentum dissipation. 


The extremal hyperscaling-violating background is given by 


ds^ = —D dt^ + B dr'^ -|- C [dxY = r ^ 


dt^ .... 

'TW + ^072" + 


{dxY" 


At ~ Qr 




( 6 . 1 ) 

Here d is for the number of space dimensions of the dual boundary held theory, ^ is 
the Lifshitz exponent, 6 is the hyperscaling-violation exponent, ( is the conduction 


exponent, |13, ^ 321 that controls the gauge held and Q is proportional to the UV 


charge density. Bq is a positive constant that can be absorbed into a redehnition 
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Figure 10: Conductivity for J = 10^. Continues lines show the Drude fitting. Gray 
dashed line is a htting with exponent —1/2 


of the generalization of the AdS scale (set to 1 here) and the bonndary space-time 
coordinates. 

It is a solntion to the EMD action, [p| , 


^A = 




( 6 . 2 ) 


where the backgronnd electric held generated by the boundary charge density is 

Qy/DB 




ZC-2 ’ 


(6.3) 


In generic scaling solutions, the bulk gauge coupling constant and scalar potential 


Z((/) ~ r'^ , 1/(0) ~ r-P 


(6.4) 


are r dependent as they depend on the running scalar. In the case where in the IR 
Z ~ e'’"^ and V ~ the dilaton 0 behaves in the IR as 0 ~ alogr therefore 
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K = a'y, p = a6. Then, the previous equation determines the conduction exponent, 

nn 

( = d — 6 — K-{ —(6.5) 
a 


g, 1^, 0 to be 


Note that this relation is only valid if the U(l) symmetry is unbroken. 

The exponents 6 and ( track the presence of the violation of naive scaling in the 
classical solution. When 6^0 the metric violates hyperscaling. When ( d, the 
gauge held prohle violates naive scaling. 

When 0 = (C = 0 there is a genuine Lifshitz scaling symmetry in the theory 
and the various observables have canonical scaling. For example the energy scales 
with mass-dimension spatial momenta with mass-dimension one, and the charge 
density p with mass dimension d while the current density has dimension z + d — 1. 
Moreover, all observables of the dual quantum held theory have naive scaling and 
the generalization of the AdS scale never appears in QFT formulae. 

This is to be contrasted with the case where at least one of the exponents 6*, is 
non-zero. In that case the following are true: 


1. There is always a scalar held that runs in the solution, 0 = f [j) and it is 
typically responsible for the nontrivial values of the exponents 6, (. 

2. Physical observables are still scaling but they violate hyperscaling and/or naive 
scaling^^ the scale £ appears in physical observables, to correct for the unusual 
dimensions. For example the entropy is 

Sr^£-^T^ , de = d-e, (6.6) 

with de appearing as the ehective spatial dimensionality. In this way S has 
physical dimension d. 


In some string theory examples, the hyperscaling violation scale £ can be inter¬ 
preted as a Kaluza-Klein scale, II , namely the radius of an internal dimension that 
has been compactihed. In this case the higher-dimensional theory is a theory that 
preserves hyperscaling. 

In the generic case, the hypescaling-violation scale £ may be distinct from other 
UV scales that drive the flow of the theory. In the simplest possible gravitational 
example with a conventional relativistic (ie AdS) fixed point and a hyperscaling 
violating solution in the IR, we must have a potential in (|6.2|) with two exponentials 
and a gauge coupling function Z with one exponential. In that case the dynamical 


define I in general by substituting r —> rji in the metric (13). 

^^Note the here we take a more general definition for hyperscaling violation compared to the 
condensed matter literature: any physical observable that scales with an anomalous exponent, 
violates hyperscaling according to our definition. 
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UV scale Kjjv can be mapped uniquely to the IR hyperscaling violating scale 1. In 
all hyperscaling-violating solutions we will therefore denote ^ by 

In the near-extremal case there is also a blackness factor /(r) in the metric 


ds^ = r d 


f{r)dt^ 

lyt 2 ^ 


+ Bq 


dr"^ 


fir) 


+ 




/(r) = 1 - (6.7) 


as well as an appropriate modihcation of the gauge held. 

From the Einstein equation we conclude that the only way for the gauge held 
to support^"^ the gravitational IR background is to have -|- C = 0. If de -|- C < 0 
the charge will not backreact to leading order and the IR solution will not depend 
on the presence of this charge. Finally, the case dg -|- C > 0 is not allowed as it 
is incompatible with the equations of motion (the gauge held contribution is larger 
than gravity). We therefore have the general inequality 


dg -\- C, = d -\- C, — 9 < 0 


( 6 . 8 ) 


Note that this inequality implies for EMD solutions that in the presence of a non¬ 
trivial charge density only the AdS 2 solution is truly scale invariant. All other solu¬ 
tions have hyperscaling violation. 

The huctuation equation for a time-dependent gauge held 6Ai = ai{r)e^^^ is 


(after solving for the associated metric perturbation, Sgu, HIT 


z-^c^ 



B , 


Q^B 

ZCd 


ai = 0 


(6.9) 


The last term is proportional to the charge density, as the charge Q enters in the 
invariant combination, —^- 3 -, which is the charge density at any given r. It is re- 

C{r)^ 

sponsible for the presence of a non-trivial Drude weight in the absence of momentum 
dissipation, |^, . 

In the near extremal scaling background (|6.1|) the huctuations obey the equa¬ 


tion 


15 


a" + 


/ 


n) + 


5, 


0 ( 

P 


f 




a* = 0, 


( 6 . 10 ) 


We will absorb the positive constant Bq into a redehnition of ca and Q, therefore from 
now on we set Bq = 1. We will also set the temperature to zero (that sets / = 1 in 


^'*This means that the contribution of the gauge field is of the same order as the gravitational 
contribution to the bulk equations of motion. 

^®It was shown in 0 that near extremality, the equation can rewritten in-terms of the scaling 
variables r/r/j and ui/T showing that the conductivity, up to an overall scale is a function of uj/T. 
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This equation is of the form 

(e"'l’'la')' + - G(r)) a = 0 


( 6 . 11 ) 


where a prime stands for a derivative with respect to r. We define a new variable u 
and a new function 'ip so that 

du 


dr 


= e 


B 


a = e-5(^+^) 


tp 


( 6 . 12 ) 


The equation ( |6.11 ) now becomes a Schrodinger-like equation 

—'ip + V-ip = uj‘^ip (6.13) 

with 

V = ]^{A + B) + -^{A + Bf+ G, (6.14) 

where a dot is a derivative with respect to u. For all gapless solutions described here 
the IR is either at r —)■ cxo or r —)■ 0 and the Schrodinger coordinate u ~ it will 
behave always as u —)■ oo. 

We now turn to our concrete problem in (|6.10|) , where the Schrodinger potential 
can be calculated asymptotically to be 
,2 _ 1 ,02 


V = 


1/2 = 


z/ 


+ 








zu 


(20 + d^-d(K + 0 + 2 )) (20 + d(de + 2 z-2-k)) 1 


Ad?z‘^ 


+ 


(6.15) 


(6.16) 


The leading form of the potential in the IR depends on the value of the charge density 
Q and the exponents. More precisely: 


I dg + C < 0 or Q = 0. In both of these cases the leading form of the potential 
will be 

— 1 

R = (6.17) 

and is scale invariant. The small u behavior of the AC conductivity is, |^, 

cT(a;) ~ cj™ , m = 2|z/| — 1 (6.18) 

II de + C = 0. When this condition is satisfied the second term in the potential 
in (|6.15|) has the same u~‘^ behavior as the hrst term and in this case 

02 

<j(u:) ~ a;"* , m = 2|z/| — 1 , z/2 = z/2 -|—- (6.19) 

z^ 

Note that in this case, the gauge held is contributing non-trivially to the IR 
equations of motion. Also the term appearing in the exponent is not the 
UV charge density (that is a free parameter) but the IR charge density which 
in this case is hxed completely and is a function of z^ 0, [^. 
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Ill The remaining case he + C > 0 is incompatible with (| 6 . 8 |) imposed by the 
equations of motion as we argued previously. 


In view of this classihcation we may now specialize 


to two different cases: 


The charged IR solution of []^ . In this case the IR background is conformal 

to Lifshitz with hyperscaling violation. The gauge held has an action 
but 7 is determined in this case by the exponents z and 6, |^. The dilaton is 
given by 

rsj 'fdS 


6 ^ = 


26 

with 

( 6 . 20 ) 



( 6 . 21 ) 


and the charge of the solution is given by 

2{z-l) 


Q = 


diQ — 1 z 


( 6 . 22 ) 


In this case the Gubser criterion reads 


z - do 


> 0 , 


dg — 1 T ^ 


d{z — 1) — 6 ’ d{z — 1) — 6 

and the thermodynamic stability condition 


> 0 , 


(^- 1 ) 


d{z — 1) — 6 


> 0 , 


> 0 . 


(6.23) 


(6.24) 


{d-l){z-l)-e 

when satished it is correlated to the existence of a gapless spectrum at zero 


and hnite temperature, \12 
For these solutions of EMD 


K = a'j = 2d — 


29{d-l) 

d 


m = 


2 dg 

3-- + — 

2; 2; 


- 1 . 


(6.25) 


In this case, both terms in the potential ( |6.15| ) scale as the IR charge density 
Q is hxed from the equations of motion and is a function of the exponents z, 6 
and we are in case II above. 

A few special cases of (|6.25|) deserve to be mentioned. The AdS 2 case is obtained 
when 2 ; —)■ 00 . In that case m = 2. The hyperscaling violating, semilocal 
geometries are obtained from the limit —)■ cxo, 0 —)■ — 00 , with ® = —7 held 
hxed. For these geometries. 


m = |3 + ? 7 | — 1 = 2 + 7 > 0 


(6.26) 
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Figure 11: Contour plots to illustrate the region in the parameter space where the ex¬ 
ponent m takes negative values for the single charged model. Left; Conductivity in the 
charged case for d = 2. Right: Conductivity in the charged regime for d = 3. The allowed 
values for the parameters are bounded by the gray mesh. The negative values for m are 
outside the permitted region. 


In Fig. ^ we plot the exponent m as a function of z, 9. We observe that the 
region where m takes negative values in the 9 — z plane, is always outside the 
region bounded by the Gubser criterion and the thermodynamic stability, for 
both d = 2,3. This was already observed in m which also gave a plot of the 
allowed values of m. It is therefore not possible to have a negative exponent m 
in this case. 


• The only other possibility that remains to be studied is when the second term 
in (|6.15|) is subleading to the hrst term. This happens when the relevant gauge 
held is subleading in the IR equations of motions and does not backreact to 
the IR geometry to leading order. To do this we will need at least two gauge 
helds in the holographic theory, one to seed a general scaling geometry and the 
other to be sub leading (and therefore treated as a probe in the IR geometry)^®. 
We introduce therefore the gauge helds A^,V^. will be nontrivial and will 
support the hyperscaling violating solution as in the previous cases. will have 
a smaller coupling in the IR and will be therefore subleading in the equations 
of motion. 


^®Clearly this generalizes to several gauge fields, with similar quantitative conclusions. However 
there could be other possibilities for seeding the leading solution that we do not address in this 
paper. 
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The gauge field action is now 


= -\j [Za{^)FI + Zv{(j>)F^] (6.27) 

and the scaling of the conductivities has been extensively analyzed in appendix 
We report the results here and we distinguish two cases. 


1. One of the gauge fields (say A) supports the IR geometry while the other 
(say V) is subleading in the IR. This was analyzed in |h.l| . 

In this case the conductivity cta is as the previous example, ( |6.25|) The 
conductivity for V, cxy scales differently. We define 


in particular we parametrize k in the following way 


(6.28) 


K = K + 6k. 


The charge exponent for V is 

- ‘26 

C = de-R + — , (6.29) 

The exponent m in the conductivity ay now reads 


m = 


d{d-2) 1 .^ ^ . 

-;- {d + z — 2 — k) 

dz z 


1 = 


z + ( — 2 


1, (6.30) 


9 and z must obey the Gubser and thermodynamic stability bounds ( |6.23|) , 
(EH as well as dg + ( < 0. 


Note that the exponent m, written in terms of the two hyperscaling vio¬ 
lation exponents 6, does not depend on 6. Therefore it is sensitive only 
to hyperscaling violation originating in the charge sector. 


A few special cases of (|6.3CI|) deserve to be mentioned. The AdS 2 case is 
obtained when z ^ oo and m = 0 in this case. The hyperscaling violating, 
semilocal geometries are obtained from the limit 2; —)■ 00, 6 —)■ —00, with 
- = —7] held fixed. For these geometries. 


d — 2 , , , 

m = —;— T] , T] > 0 , d>2 (6.31) 

d 

and m is always positive. 

In figure |T^ we show the values of the exponent m for d = 2 in a charged 
background {Qa, Qv 7^ 0) where the backreaction of Qv can be neglected. 
In general m depends on z, 9, R. 
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Figure 12: Contour plots to illustrate the region in the parameter space where the expo¬ 
nents m takes negative values for the V gauge field. The plot was made for d = 2. The 
negative values for m are inside the permitted region. We have set |(5 k| =3 


The plot represents m for \R — k\ = 3. It is always possible to obtain a 
region with negative values of m, either with R — k > 0 and z > 0 or 
R — K, < 0 and z < 0 (see appendix 0 . Note that from Eq. ( |6.30|) m = 0 
for R = 0 which is the case when the bulk gauge coupling constant flows 
to a constant value in the IR. 

We conclude that in order to have m < 0 in d = 2 two conditions must 
be satished: 

(a) The charge density in question must not backreact on the metric. 

(b) k > 0, namely the bulk gauge coupling must be IR free for 2 ; > 0 and 
UV-free for z < 0. 

The inverse correlation to the sign of 2 ; is also correlated to the relative 
flow between energy and momentum, controlled by — ~ ,, 2 ( 1 - 2 ) taking 
into account the constraints this ratio for z > 1 vanishes in the IR while 
for z < 1 it diverges in the IR. 

2. A different case involves the two gauge helds being equally important in 
the IR, but that their coupling constants, Za and Zy, become different as 
one flows to the UV. In that cases Za — Zy is much smaller than Za + Zy 
in the IR. This cases was analyzed in E3- What was found is that a a and 
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cry contain at low frequency, a scaling term of the type shown in ( |6.25|) 
and another as in (|6.3CI|) . There is also a cross-conductivity cr^y that is 
partly negative and contains the difference of the two terms above. 


7. Outlook 


We have analyzed the AC conductivity in a holographic model of a strange metal 
proposed in [0. We found some interesting scaling limits of the conductivity that 
seem to generalize to general critical holographic systems. We have also found that 
this scaling regime is controlled by the uncharged system and seems to be generically 
independent on the mechanisms of momentum dissipation. 

The detailed results are as follows: 


1. We dehne a parameter q that distinguishes the Drude regime (g S> 1) where 
the DC conductivity is dominated by dissipation and the Charge Conjugation 
symmetric regime {q <C 1) where dissipation is negligible. 


2 . The characteristic temperature scale that controls the AC conductivity is an 
effective temperature Tg// (with an associated scaling effective temperature fg// 
dehned analogously). This is a characteristic effect for the system in question 
and more generally for systems where charge is described by a D-brane system. 
It is not expected in more general cases. 


The effective temperature found is distinct from the system temperature T and 
is due to the existence of a hierarchy of interactions. Similar effects have been 
observed in electronic systems and in holography, [^-[^ . 


3. We have dehned a generalized relaxation time r by the IR expansion of the AC 
conductivity, 

"ze aoc {^ + i T ^ + 0{ijj‘^)) (7.1) 

In the presence of the Drude peak, this is the conventional dehnition of an 
associated relaxation time. When there is no Drude peak present, r is still 
well-dehned, although in that case the interpretation as a relaxation time is 
not clear. 

We gave an analytical formula for r in ( [4.14|) . It takes a simple form for large 
and small values of the scaling temperature variable t. In the regime I (see hg. 
H) we obtain 

r ~ \/T (7.2) 

while in the regime II (with t typically large) it is set by the inverse of the 
temperature 

r ~ 1 (7.3) 
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4. In the Drude regime {q 3> 1) where the dominant mechanism for the condnc- 
tivity is dissipative, there is a clear Drnde peak as seen for example in hgnres 


In the Charge Conjngation Symmetric regime it is also clear from the same 
hgnres that no Drnde peak is to be seen in the IR of the AC condnctivity. 

5. At zero charge density (Charge Conjngation Symmetric regime) there is a scal¬ 
ing tail for the AC condnctivity that behaves as 


cr 


u 


"eff 


71 


Arg{a) ~ - 
0 


(7.4) 


For hnite charge density this tail survives not only in the Charge Conjugation 
Symmetric regime but also in part of the Drude regime, as seen in the various 
plots of hgure ^ as well as the ones of hgure ||. The qualitative reason for this is 
that in the presence of the Drude peak, its tails falls off as ^ and this is faster 
than u~3. Therefore, the scaling tail will eventually win over the Drude peak 
for uj > ujo and the only condition that it is visible is that the UV structure of 
the theory kicks-in at uuv 3> ojq. 

6 . This scaling tail of the AC conductivity generalizes to more general scaling 
holographic geometries, as previously described in O- The equation that 


determines the conductivity is given in (|6.9|) and the equivalent Schrodinger 
problem has a potential of the form Vefj = Vi p^V 2 where p is the IR charge 
density and is proportional to the UV charge density. 

In particular for a metric with Lifshitz exponent z, hyperscaling violation ex¬ 
ponent 9 and conduction exponent C, with d spatial boundary dimensions, we 
hnd that in general 

TT m 


cr 


with 


UJ 


m = 


Arg{a) ~ 


z + C — 2 


- 1 , 


(7.6) 


(7.6) 


There are several constraints in the parameters of this formula that are detailed 
in section 1^. 


7. There are some special cases of ([7.6|) that deserve mentioning. For an AdS2 
IR geometry the exponent can be obtained by an z —>■ cx) limit in ( [f.6|) giving 
m = 0. 

For hyperscaling violating semilocal geometries we must take 6 ^ oo, z ^ oo 
with - = —g hxed and obtain 


d- 2 


m = 


-g + l 


1 = 


d-2 


d 


g 


(7.7) 
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Finally, for the gauge field conformal case we obtain m = 0 when d = 2. 

8 . We find that for two spatial dimensions, negative values for the exponent m 
are correlated with the sign of the Lifshitz exponent z and the strength of the 
gauge field interaction in the bulk. Parametrizing the IR asymptotics of gauge 
coupling function Z as 

Z r'^ , r —)■ cxo (7.8) 

in conformal coordinates, we obtain that zk > 0 for negative values of m to be 
possible. 


9. In the special case where the associated gauge field seeds the IR scaling geom¬ 
etry, K is fixed as a function of z, 6 and the exponent m takes the value 


m = 


3z-2 +d-e 


- 1 


(7.9) 


and is always positive. In the case where the IR geometry is AdS2 can be 
obtained for z ^ oo limit of (|7.9|) giving m = 2. For hyperscaling violating 
semilocal geometries, we must take, 9 ^ oo, z ^ oo with ® = —rj hxed. In 
this case we obtain 

m = \3 + T]\ — 1 = 2 + 7] (7.10) 


10. An important issue is whether the scaling of the AC conductivity described 
above for the general scaling geometries is controlled by the dynamics of the 
charge density, or it is decided by the neutral system. 

In the example we have analyzed the effective Schrodinger potentials that con¬ 
trol the calculation of the conductivity have two parts. One that is independent 
of charge density and one that is proportional to the square of the charge den¬ 
sity. In the generic case it is the hrst that controls the UV scaling of the AC 
conductivity described above. Only if the charge density is supporting the IR 
geometry, then the second part is of the same order as the hrst and it is the 
sum that controls the scaling of the AC conductivity. This special case is also 
the only one we found where the exponent m in (|7.6|) is always positive. In 
other cases it can also be negative, but unitarity implies always that m > — 1. 


These hndings suggest that there is a generic source of scaling tails in the AC 
conductivity in holographic systems. Moreover, in generic systems this scaling is ex¬ 
pected to be independent of the mechanism of momentum dissipation. On the other 
hand, whether it is visible it depends on the details of the momentum dissipation 
mechanism. 

An important part of this story is the dependence of the T=0 result on scales. 
For this we must consider a hierarchy of cases as a function of the properties of the 
system. It also matters whether the charge density is interacting strongly with the 
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gravitational system or is a probe whose backreaction on the system can be neglected. 
We will discuss all these cases in turn 


In the truly scale invariant hxed points (C = 6* = 0) with a general Lifshitz 
scaling exponent, z, the AC conductivity will scale with u with the canonical 
dimension, 

a ~ (7.11) 

for any z and no intrinsic scale enters in this relation. Note that in the EMD 
case, apart from the AdSd+2 solution (z=l), there are no other such scale- 
invariant solutions. To obtain general scale invariant Lifshitz solutions, one 


should use the more general EMD action, |^] where the U(l) symmetry is 
generically broken. 

In the presence of hyperscaling violation 6^ 7^ 0, and/or C 7 ^ Oj the IR theory 
contains a nontrivial scale A/r that appears in many relations and is responsible 
for the change of scaling exponents and the violation of hyperscaling, 0- In 
some cases which can become regular by embedding the dual bulk theory to a 
higher-dimensional relativistic scale invariant theory such a scale is the radius 
of the extra dimensions. This scale appears in scaling relations in order to 
“correct” the naive scaling dimension. For example 


(j{uj) ~ A 


d—2—m 

IR 


u' 


(7.12) 


• In the cases where the charge density acts as a probe in the IR geometry (and 
its effects are therefore subleading), it introduces a new IR scale that is set by 
the associated charge density. However we have seen that this new scale does 
not enter to leading order in the scaling relation for the associated conductivity. 
We are therefore again in the situation described in the previous cases. 


It is expected that hyperscaling violating geometries should have a resolution 
of the their naked singularities via an AdS 2 or a stringy geometry in the IR, 
and they will be only intermediate scaling regimes, ||^. In such a case A/^j is 
the associated scale to this mid-IR regime. 


To try and extend qualitatively our general scaling analysis of section ^ we must 
reason as follows: There are two further modihcations to the setup above so that 
this result is embedded in a complete RG flow. The hrst step is to match these IR 
asymptotics to an AdS UV geometry. The second it to turn on hnite temperature. 
We will now discuss such modihcations in turn. 

The non-trivial RG how at T = 0 introduces in the simplest of cases a transition 
scale Auv In the hyperscaling violating geometries this is related to the characteristic 
scale of the IR geometry. In scaling geometries it is related to the UV operator that 
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drives the flow. Quite generically this can be the charge density. Our general AC 
scaling found above is valid in the region uj <C Auv- 

In the presence of a finite temperature, T, it was shown in [|^ that the AC 
conductivity is a function of oj/T. There are two ratios that control the behavior of 
the AC conductivity: T/Ajjv and uj/T. In the hyperscaling violating case, Air also 
enters. However, in the limit Ajjv oo. Air enters trivially as in ( [7.12|) . We have 
the following regimes 

• S> 1. The theory is in the UV scaling region and the AC conductivity is 
characterized by the UV fixed point alone. 

• 1. In this case, there are three regimes for the AC conductivity as a 
function of frequency. 

— a; -C T This is the IR regime controlled by the properties of the black- 
hole horizon. In this regime the mechanisms of momentum dissipation and 
hydrodynamics are controlling the behavior, |^. If the drag dissipation 
is present, it determines the Drude peak (that maybe pronounced or less 
prominent depending on the strength) and this is the main characteristic 
of this region. 

— T -C a; -C Auv In this regime the scaling behavior found earlier, at T = 0 
is prominent. As the power evolution of the AC conductivity is always 
~ oo"^ with m > —1, it will always dominate over the Drude l/uj decay for 
large enough frequency. The only possibility where this will not be visible 
is if the cross over scale is of the order of or larger than Auv- 

— o; 3> Auv- In this regime the AC conductivity is determined by the UV 
fixed point. 

The type of AC conductivity scaling described here is very reminiscent of the 
one seen in the cuprates in in the intermediate region T -C a; Auv- Indeed 
the scaling exponent found there is m = — | which suggests in view of our discussion 
the presence of hyperscaling violation. 

There are several problems that remain open in view of our results 

1. We should back up our scaling analysis and the description of different regimes 
above with a complete calculation of the AC conductivity from numerical com¬ 
plete RG flows that connect AdS UV fixed to either IR or intermediate geome¬ 
tries with general scaling exponents {z, 9, (). 

2. The analysis of scaling described here was based on geometries that preserve 
the U(l) symmetry. This analysis needs to be extended to the general class of 
U(l)-breaking extremal scaling geometries found in |^. Part of this analysis 
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was already done in but it should be extended to more general cases. We 
expect that our results will extend smoothly to these more general solutions. 
The aim is a universal formula for the AC exponent m that controls the scaling 
of the AC conductivity. 

3. Although the scaling tails described above are expected to be independent 
of the mechanism of momentum dissipation, it is interesting to study their 
interaction with momentum dissipation and different types of low-frequency 
behavior. In particular it would be interesting to study this in geometries with 
helical symmetry as the general class presented in ^ . 


4. The scaling exponent of the AC conductivity is expected to depend on the 
charge interactions captured by the bulk (self)-interactions of the U(l) gauge 
bosons. Understanding how is an interesting question. 


5. The eventual comparison with cases realized experimentally as in 
interest. 


is of direct 


We are planning to address some of these questions in the near future. 
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APPENDIX 


A. Effective Lagrangian and Equation of Motion for Gauge 
fluctuations 


In this appendix we present the technical details of the derivation of the equation 
that determines holographically the two-point correlator for the current and the 
associated AC conductivity. 

We substitute for the gauge field the background value as well as a linearized 
perturbation A = A(o) + ea. e 1 is a small parameter that we use to derive the 
linearized perturbation equations. 


£ = -M^/- det (7 ef) 

~ -M\/- det 7 (^1 + |Tr [7"V] + ^Tr [y’V] ^ - ^Tr ^ ^ (A.l) 

where we have defined 

7 = ^ + ^(0) , / = da . (A.2) 

We can split the matrix 7“^ into its symmetric and antisymmetric parts 7“^ = 
respectively. The symmetric part of 7“^ is as usual interpreted as the inverse open¬ 
string metric. This induced world-volume metric has a horizon at the surface u = m* 
instead of u = Uq, and it is given by 


s — g — F{Q)g ^F{Qy (A.3) 

The quadratic part Te// in e of ( |A.1|) will determine the dynamics of the fluctu¬ 
ation a. It can be rewritten in the form of a Maxwell theory with a Chern-Simons 
contribution as follows 


Jeff — —MV— det s + 


8\/— det s 


M^^^^fMNfpqQR] , (A.4) 




where the effective coupling is 

= V— det s/ \/— det 7 


9 I 


txyzu _ _ _ 1 

^ — ^txyzu — -L 


nMNnPQ 
^MNPQR^ ^ 


(A.5) 

(A.6) 


Notice that now we are raising indices with the open string metric s. In particular 
the induced line element is 

ds^ = + ‘2s+-dx~^dx~ + 2s+udx^du + 2s+idx~^dx^ + s _ {dx~Y + 

2s-idx~dx^ + Sijdx^dx^ + SuudM + 2suidudx ^, (A.7) 


The components of the world-volume metric above are given in appendix ^ in terms 
of the background metric. Finally, the equation of motion for the perturbation is 


9„ (= 0, (A,8) 
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B. The Open String Metric 


In this appendix we give the explicit expression of the open string metric in terms 
of the background metric and the other background fields. The charge density and 
the background electromagnetic held modify the induced geometry on the probe 
D-branes. Using eq. (p.4|) we can write the non-vanishing components of the world- 
volume metric as 


I UP h'l 
S ++ — g ++ H - 1 - 


h' h' 


9i' 


9u 


S+- = 9+- 




Eyh’y + 

9ii 


h[h'. E,Ejg__ 
Sij — 9ij H-1" 


9u 






+ ■ 



Quu 



+ 

— 


Quu 






h'+K 


S-i = 


Quu 

i'K 




— (^+-h'_ - 9—h'+) 


Quu “ 1 “ 


+ , 9 -( fc '+)" + 9 ++( fc -)"- 2 g +- fc '_/4 


9i' 


+ 




(B.l) 


with i = y,z and where we use through all the appendices the notation E = {Ey, E^), 
h = {hy, hz), E ■ h = Eyhy + Ezhz and E x h = Eyhz — Ezhy. Please do not confuse 
h with the blackening factor h{u) (see Eq. O). 

The components of the antisymmetric matrix are 

e+- = -^^E-h' , - g+_h’_), 


det 7 


det 7 


e+y = E. 


e+^ = E 


gu {g—guu + hP) g—h’^ 



9-y = U 


0-^ = E 


det 7 det 7 

gu {g+-guu + h'_h'^) g+-h'^ 



KK9- 

^ det 7 ’ 

_ > E'yh'yg— 

^ det 7 ’ 

Ez 9 +-h'yK 
det 7 ’ 

E y9+-h'zh 'y 
det 7 


= 


E\^9iih- , 2 (9++h- 9+-h'+ 

— -+ g.. I ^-^- 


det 7 


\ det 7 det 7 


^ EyEzg—h', ^ Elg__h'y ^ G+_guh'y 
det 7 det 7 det 7 

Qzu _ EyEzg — ^ Eyg — hy ^ ( j-i — gah^; 

det 7 det 7 det 7 

E xh' {g+-h' - g—h'^) 


Qyz = 


det 7 


(B.2) 
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Finally having the components for s and 6 we can compute the vector Q (|A.6|) , which 
reads 

Q = —E X h'{g+_dx~^ + g _dx”) + guh'_E x dx. (B.3) 


C. Change of Coordinates 


In this appendix we show the f's functions appearing in the change of coordinates 
used to diagonalize the open string metric (|3.4|) . For completeness we write again 
the transformation 


dx"*" = dr + f^(u)dX + /+ (M)dM , 
dx_ = dX" + f;{u)dT + f~ {u)du , 
dy = dY + fr{u)dT + fl{u)du , 

dz = dZ + /;(w)dr + f:{u)du , (C.l) 

with the radial coordinate rescaled as m —)■ The functions f's read 


fu 

n 

fr- 

fu 


b\h-l) 


Ui,E ■ h' 

h {h^E^ + gl^ - b‘^E^ 


u^E-h 

+ 1)—7—r-\ 

^ h (b^E^ + gl \ - b^E^ 

h'_ 

K 

h' - glih - 1)) + 2ffgi{h + l)h\ 


262 gg (262(ft - l)ftV - (h + l)h'_) 

^ M* {{h + l)h'_ - 2b‘^{h - l)h'^) 

2{b^{h-l)E^ + glh^ 


(C.2) 

(C.3) 

(C.4) 

(C.5) 

(C.6) 


. 2u^h - l)(?i+)^ - (h'_)^ (^ib^\E\^ - gf-(h - 1)) - Ab^gf^(h + l)h'_h'jd^ (b^\eijEj\(l - h)E x h' + 9 iihh'^'j 

b'^g?. ((h + l)h'_ - 2b'^{h - (h (b'^ + 4 (^E x h'Y u3(h - 1)^ + g?. fA(h')'2u^h + l)^ - 62 |£;| 2 ^ 
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After using this transformations the new components of the metric can be written 
in term of the old one as follows 


Suu — ‘^fu {^+-fu + ^-yfu + ^mzfu) + 5- {fu Y + “^fYi^+yfu + ^+zfu + Su+M*) + 

S++{fYY + ‘^fui^yzfu + SuzU^) + ‘2SuyUjy + Syyif^y + + SuuUl, 

(C.8) 

■Srr = 2/^ {s-yfr + ^-zfr + -5+-) + ^-(/,- + 2f^{Syzf^ + S+j^) + 

^yyifrY + ‘^S+zfr + SzzifrY + -5++ ) (C-9) 

S = 2/l^(S + y/^ + S+zfl + S + _) + S + +(/l^)^ + 2f'^{Syzfm + S-y) + Syy{f^Y ^ 

+‘2‘S-zfl + SzziflY + ■S-; (C.IO) 

and 


^yy ~ ^yy > 

5 ^ 2 : — • 


(C.ll) 

(C.12) 

(C.13) 


D. Fluctuations 

This appendix provides the details on the computation of the equation of motion for 
the fluctuations. The solutions of this equations are necessary to obtain the retarded 
correlator among two charged currents. 

In order to do so, we will set Ez = 0 and we will introduce the fluctuations 
ttM^u, x~, T, Y) into the Eqs. (|3.5|) . With this dependence in coordinates the system 
can be decompose in terms of longitudinal and transverse modes respect to the y 
direction. We can also fix the gauge a„ = 0. Thereafter, we Fourier transform them 
as follow 

f{u, X-, T, Y) ^ ^ (D.i) 

where (3 = After some tedious computation we get a set of three differential 

equations plus one constraint in the longitudinal sector {t,—,Y) and one differential 
equation in the transverse one (z—direction). The constraint read as 


— ks'^^a'y{u) + oj (s'^’'a).(M) — a'_{u)) = 0 , 


(D.2) 


have identified x'^ as the boundary field theory time, however the new temporal coordinate 


and the former time are related at the boundary through the transformations (3.6) as dz^ = dr — 
3lPdX~. Therefore, we evolve in the Fourier transformation along the direction of a;+ = t — [ 3 h^ X~ . 
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the longitudinal equations are 


ci'Au) + a!_{u^ 


2 / ggUUg 


)' 






iMQ':, 


—I3roj^SuuS^^o^riu) + l3b^kijJSuuS^^CLy{u) - ^- S Suu {ujay{u) + kar{u)) = 0 , 

£/5 V-s 

(D.3) 


a"(M) + a'^{u)^^ 


■2 


"V—ss““s'^’' 


— ar{u)Suu{b‘^'l‘^^‘^S + + 


-f3b‘^uj'^SuuS a-{u) — kojSuuS^^ay{u) — ^jl^ -SuuSTT {ka-{u) — f3b'^ujay{u)) = 0 , 




:uu ~yy 


ni ^ . n ^ (^5 

ay[u) + ay(M)— 32- 


)' 


(D.4) 


/ ~„y ay{u)u^~Suu{li^b^~s + 5 "") + 

g^ y/—ss^'^syy 

b^f3kuSuuS a-{u) — kuSuuS^^ariu) + gn„g^^(a_ {u) + (3b‘^ar{u)) = 0 , 

95 V-s 

(D.5) 


and in the transverse sector we obtain 

[(3'^b'^s + s’"^) + az{u) = 0 , (D. 6 ) 

notice that the Chern-Simons term vanishes in the last sector, because if = 0 the 
only non-vanishing component of Q points in the z direction 


Q EyQiih_diZ. 

(D.7) 

Now it is possible to eliminate one of the equations in the longitudinal sector 
using the constraint (U. 2 ) and the gauge invariant helds, ending with the same 
number of equations as propagating degrees of freedom. To do so, let us redehne the 
helds in term of its gauge invariant electric helds. 

S_{u) = — [a_{u) + (]b‘^ar{u)) , 

(D, 8 ) 

k 

Sy{u) = ay{u) H— ar{u ), 

(D.9) 

Sz{u) = ttz 

(D.IO) 

now using the following combination for the equations 


-(M- + ^b^M^) , My + 

Ul 

(D.ll) 


and solving for a'y^u), a'_{u), a'^{u) from Eqs. (E3). (E3) and O we obtain the 
two linearly independent equations for the longitudinal sector 
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(b'^'yojclc + c_(kc^ — uc'^)) , , ^ , 9 _ 

^-('“) + 7^r + d_)8-{u) 


uj{b'^'yc~ — c'^) + key 


{b^-fdy + dl) Sy{u) = 0 (D.12) 


c"r ^ , -c") + kclcy) \ c / N , 

(“) “I 77T8Z^ ^kWT^, + {~ + ^y] 4(“) + 


uj{b‘^'yc~ — o’") + 


a; 


Zee {cl-cl) 


{kudl — k{kdy + udy) + oj’^d'l) 


~^ £-{u) + '" ^~^' s.{u) = o (d.i3) 

oj{b^'yc~ — c"^) + key b^'yu^ 

where the unknown functions c's and d's are the coefficients of the equations of 
motion and constraint, read it from Eqs. (p.2|) - (p.6|) as follow 

a”{u) + cia'j{u) + d{aj{u) = 0 , (D.14) 

cPa'j{u) = 0, (D.15) 

with i,j taking values —,T,y. Notice that the equation for Sz is no written because 
it remains the same as Eq. (p.6|) 


E. Sources and Perturbative Solution 


In what follows we will explain the procedure to obtain the perturbative solution dis¬ 
cussed in the subsection M This solution was useful to obtain the DC conductivity 
and the relaxation time Eqs. ( p.l3D , ( [4.14D . 

First of all we plug in the ansatz ([4.61) into the Eq. ( [4.3|) , so the system can be 
rewritten as follows 


du {aoiu)a{u}duA^^\u)) ^ 5*’ ^\u) , 


(E.l) 


where the index i takes values i = 0,1,2. The particular form of the sources up to 
second order in frequencies are 


2N' 


iS® = idu ((1 + u)a{u)A^^'^ + (—1 u^)a{u) log (1 — u)X^ 


( 0 ) 


(E.2) 


= du ((1 «)«(«) log(l - 4- -{-l + u^)a{u) log(l -uf + 


i{l + u)a{u)X'P + i(—1 -4 X)a{u) log (1 — u)XX 




eff 

T 


with the zero-th order source Si = 0, and where we have dehned 


p{u) = -2g_ 




bHMu 


P^- 


p{u)u 
i-l + U^) ’ 

(E.3) 


(E.4) 
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The general solution for eq. (^) can be written as 


4‘>(t.) = C® + [ di. ' (c® + [duSt»(u)) , (E.5) 

J ao{u)a{u) \ J J 

the integration constants are hxed demanding regularity for the solution at the 
horizon and is directly identihed as the non-normalizable mode. 

If we start solving order by order, we get for the zero-th order solution 

Af > = C|”> = a,(0), (E.6) 

(72°^ must vanishes in order to have a regular solution at the horizon. At the next or¬ 
der the source S doesn’t vanishes, in consequence, regularity at the horizon demands 
that 




■Kl^Tefful 

Finally the hrst order solution can be written as 


A^p{u) = 


log (1 — m)- — log 


_a(l) 


a{u) 


1 — u 


“d^Ahiog 


1 + u J Jq a{xy 


If we repeat the same analysis at second order, regularity demands that 


Cf ^ = a,(0) 


log(2) - 


"dx^log 


1 — X 


^27rTe//M* ' 4 Mul Jq a{xy yl-fa: 

and in consequence the solution at second order is 
1 


(E.7) 


1 — X 
1 -I- X 

(E.8) 


(E.9) 


XPin) = 


u — 1 


iAj\u) + i{u — 1) log (1 — u)XJ^\u) + 02 ( 0 ) log (1 — u) + 
27r^6^az(0)£^Mop(M) log (1 — v?) 


yA/" (1 -|- u) a{u) 

2'K%'^aziO)tTjjj^uo A 
(1 + u) a(u) 


(1 -|- u) a{u) 
dxp'{u) log (1 — x^) j . 


+ (E.IO) 
(E.ll) 


F. Two charges hyperscaling violating geometries 


We will study the following Einstein Maxwell Dilaton system with the particularity 
of having two gauge helds in order to assess the interaction of different types of charge 
and their effect on the scaling of AC conductivities. 

The action is a simple generalization of the EMD action 


S 




^l(0) r2 Z 2 X) jp2 


(F.l) 
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The equations of motions for the background (|6.1|) read 

C" DC’ fB' a D'\ ^ 

1 V- + ~ + - — + — = n 

2^^ DJ\B C d) \Z2 Zi) C D 
Ql , Ql\ B , dC f{d-l)C' , D'\ 1 2 


Z2^ Zi) 2C'^ ^^^2C\ 2C ^ d) 2^ 


(F.2) 


where Qi and Q 2 are the charge densities of Ai and A 2 respectively. The gauge helds 
are given by 


,, Ql bd 


26 I ^ 18 

The fluctuation equations for the two gauge helds are 


(F.3) 

(F.4) 


where 0 = y + 


Z^^C-^ ^ ^ lz,a’,^ ^C^j - ^ (Qiai + Qsaa) + = 0 (F.5) 


„ 1 „ d-2 D ! , D „d-2 

^2 c ^ V 5 ^ 


Q 2 D 

Z2C^ 


{Qitti + Q 2 O 2 ) + Qj‘^cb2 — 0. (F.6) 


At this point we shall try to hnd a hyperscaling-violating metric as a solution of 
the background equations. To do so, we substitute (|6.1|) into (F.2|) and obtain 


2q2 

Q? + 2Bz - 2— - 2d{z - 1) = 0 

(Jj 

2(z - l)(d„ + z)- j = 0 

Bn fA+A). 2r-‘-‘+A + 


Zi Z 2 


9 


2de((d + 2z-l)--id + l))-a^] = 0 , 

d 


(F.7) 

(F.8) 


(F.9) 


These equations can be analyzed in two simplihed regimes depending on the ratio of 
Zi and Z 2 in the IR. 


^®We have assumed Zi ^ r'^' and the background functions correspond to the hyperscaling- 
violating Lifshitz geometry (O) 
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F.l Background charged under a single gauge field {Z 2 3> Zi) 

This configuration makes the gauge field A 2 ^t subleading in the IR region and non 
backreacting in the geometry. We will parametrize the couplings Zi and Z 2 in the 
following way 

Zi = r'^ (F.IO) 

Z2 = , (F.ll) 

where 5k will be required to be positive or negative depending whether the IR is 
at cxD or zero respectively. The lagrangian parameters can be fixed in term of the 
parameters in the metric as follow. 

(7 - 5)2 + 2 d( 5(7 - 5) + 1) „ 5cP 


z = 


Qi — 


(7-5)(7 + 5(d- 1)) 


9 = 


7 + 5{d - 1) 


dg Z — 1 


K — 2dg H —- 
d 


(F.12) 


(F.13) 


5k is not fixed, but under the constraint 5/t > 0 if the IR is at r —)■ cxd or 5/t; < 0 if 
the IR is at r —)■ 0. The Gubser criterion for this solution is 

{z — l){dg + z — 1) > 0 , {dg + z — l){dg + z) > 0 , dg{dz — 6) > 0 , (F.14) 

and the thermodynamic stability condition combined with the existence of a consis¬ 
tent near extremal black-hole is 


z{dg — z) < 0 . 


(F.15) 


Having the background geometry we can write the equations for the gauge field 
fluctuations as 


a'; + fixlxhl. B„Q,Q., 




-(X\ — 




(l2 


n (3—£— —a2 + Bou'^r'^^ ‘^a2 


_ BqQ2Qi sk 


-r 02 = 


r oi 


Now we perform the following change of coordinates 

/'xz\^B 

r ^ — 

\ w / 

^-2+Ci 

tti ^ r 2 Oj, 

which allow us to write the equations in the following Bessel-like form 

2W| !, f 2 2\ 

xa^+xa^ + [x - n J Oi = — I — 1 4/202 

Z ' ZOb' 

x‘^a 2 + xa'2 + - n^) 02 = % f—) Oi -h Q2 f—) " 02 

\ / z^ \ \zx/ ^ 


\zx 


(F.16) 

(F.17) 

(F.18) 

(F.19) 

(F.20) 

(F.21) 

(F.22) 
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where 


n 


2 

1 


n 


2 

2 


[dg — z 2)^ + 4Qi 

((5/t -|- dg — z 2)^ 
4^:2 


(F.23) 

(F.24) 


We are interested in the IR behavior of the AC conductivities. With the equa¬ 
tions written in this form and using the fact that zSk > 0, we can do perturbation 
theory expanding the helds as follows 



(F.25) 


The zero-th order functions satisfy the homogeneous Bessel equations, with infalling 
boundary conditions. The solutions are 



(x) 



-3/4g^ 




X) 


(F.26) 


where Hui (x) are Hankel functions of the hrst kind. The equations for the hrst order 
perturbations are again Bessel equations, however in this case they are inhomoge¬ 
neous 

2 (1)" I (1)' I / 2 2^ (1) 

X + xal[X — al' 

2 ( 1 )" , ( 1 )' , / 2 2 \ ( 1 ) 

X a2 + xa2 -1- (x — n2 j 

with 

F^{x) = {zx)~^ af (F.29) 

F2[x) = {zx)-^ (F.30) 


= Fi{x) (F.27) 

= F2ix ), (F.28) 


The solution for this system can be found using variation of parameters, and can 
be written as follows. 



(F.31) 


At this point we are not interested in solving exactly the integrals, but in checking 
whether this hrst order solution are in fact subleading respect to the zero-th order 
one, and the perturbative analysis is not breaking down. To do so, we will use the 
fact that the IR is at a; ^ ±cxd and the asymptotic behavior of the Hankel functions 
is of the type 

Hn{x)r^x-F\ (F.32) 
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modulo oscillating functions for the amplitude. Using this scaling, we may check 
that the dominant part in the asymptotic expansion is of the form 



_3_ 5 k 


(F.33) 


which is consistent with the perturbative analysis. 

Finally, after knowing the system can be decoupled at leading order in the IR, 
we can follow the computation of [^, to calculate the reflection coefficients of the 
analogue Shrodinger problem. Hence, the conductivities of each gauge held can be 
computed 


CTi ~ (U 


(72 ^ to ' 


m = 


n = 


2 — 3z — d, 


z 

6 k 6$ — z 2 


- 1 . 


(F.34) 

(F.35) 


As the IR geometry is supported by the charge Qi, the huctuations of ai will feel this 
charge. That explains the same result for the conductivity as in the single charged 
gauge held Eq. ( |6-25 ). On the other hand the background is neutral for the second 
gauge held. Therefore, the conductivity depends on the background parameters 
z,9 and the exponent of the gauge coupling Z 2 = and the conductivity 

coincides with the one for a single neutral gauge held. 

We have plotted in hgure |^, the Gubser bound of (|F.14|) , the thermodynamic 
constraint ([F.15|) and the negative values of m (left plot) and n (right plot) for d = 2. 
The value oi 6k has been hxed to |5«:| = 3. We observe that negative values of m 
are excluded after imposing all the constraints. However the possibility of having a 
negative n is not excluded, as we observe in the right-hand side plot. 


F.2 Both charges equally backreact the IR geometry Z 2 ~ Zi 

In this case we will parametrize the couplings as follows 

^ i (Z + Z) (F.36) 

Z2 ^ i (Z - Z) , (F.37) 

and 

Z = (F.38) 

Z = r^ (F.39) 

the parameters k and k will be chosen in such a way that Z is subleading in the IR. 
In other words, after the substitution of Zi and Z 2 into Eqs. O - (E3). we obtain 
the following set of algebraic equations 
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Figure 13: Contour plots to illustrate the region in the parameter space where the expo¬ 
nents m (left) and n (right) take negative values for the second charge not backreacting 
the background. Both plots are made for d = 2. The negative values for m are outside the 
permitted region, but there is an allowed region for n. We have set |5 k| = 3 


pns 

+ 2ez - 2— - 2d(z - 1) = 0 (F.40) 

d 

2 Bo{Ql + - 2 {Z - l){dg + Z) = - 2 Bo {QI + Ql) 

(F.41) 

Bo {Ql + Ql) - + (F.42) 

(^ 2 de{id + 2z-l)- ^{d + 1)) - = -Bq {QI - Qj) r^<^e- 2 n+-K+^ ^ 

if the following constraint on R is satisfied 

29 

Cfi = K — - — 2 dg <0 , r —)■ oo (F.43) 

(Jj 

29 

K — - — 2 de >0 , r —)■ 0, (F.44) 

(Jj 

at leading order the system becomes an algebraic system and the parameters can be 
hxed for the IR geometry. They read 


= 


(7 - ( 5 )^ + 2 d (^(7 - ^) + 1 ) 

(7-(5)(7 + 5(d- 1)) 


9 = 


6d^ 


7 + 6{d - 1) 


Qt — 


do -\- z — 1 


Kj — 2dQ T 


29 

~d 


(F.45) 

(F.46) 
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where Q‘^ = 2{Ql + Q^). 

The fluctuations in this limit obey the following equations of motion 


+ ^ ^ + 5oa;V^-2ai - {Q^ai - Q2a2) = (F.47) 

rj^ rjrt ^ 

-r'" (R - k) a[ + + (5202)^ (F.48) 

02 + -—^—^02 + Bou'^r‘^^~‘^a2 - (Qifli - <5202) = (F.49) 

r'" (it - k) 02 - + < 5202)^ , (F.50) 

The form of the equations suggest the following redefinitions 

^1 = Qi(^ 2 ~ Qio^i (F-51) 

^2 = Ql®l + Q2®2 ) 


and to combine the equations as follows Qieq 2 — Q 2 eqi and Qieqi + Q 2 &q 2 - After 
doing so, we obtain the following set of equations 

A" + = r-2+'=^C>i (F.52) 

r 

A" + + Bo A 2 = r-^+^-02 , (F.53) 

where 

Oi = (F.54) 

, „2 A 2(i^ — d(29 + it)) ((Qi — + 2 Q 1 Q 2 A 2 ) — 2BodQiQ2Q^A.2) 

Qjy^rp 

O 2 = (FAS) 

{2BodA2 {Qt -Qt)-r {26 + 2d^ - d{2e + R)) {{Qj - Ql)A'^ - 2Q,Q2A[)) 

(JjKc^ rj^ 

Notice that the factor in the right-hand side of the previous equations is 

even more subleading in the IR than the right-hand side of the analogous equations 
(Eqs. ([FA6|) and ( [FA?]) ) of the previous part. Hence, these set of equations can be 
written in terms of the Bessel equation and the conclusions will be exactly the same, 
as in the previous case. We skip this analysis and we will directly work with Oi and 
O 2 set to zero at leading order. 

In consequence, both redefined gauge fields diagonalize the system and satisfy 
the same equation, however Ai is neutral and the held A 2 has charge Qt- This fact 
implies that the normalizable modes of A^^\A^'^ will depend linearly only on its 
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own source respectively, hence, = (2)742*^^® Now inverting the 

change of variable ( [F.51|) we can write the normalizable modes of Oi and 02 as follows 


W = ^ (0i(2)4“>(a<“>,4“>) - Qj{l)4“>(aS“>,4“>)) . (F,56) 

(J 2 ) = ^ 4“’) + Q2(2)4“Vr. 4°')) . (F-57) 

which at the same time can be rewritten as 

W = 7^ (<?;(2> + of’ + ((2) - (1)) 4"’. (F.58) 

W = 7^ ((?42> + Q?(l)) 4"’ + ((2) - (1)) 4"’ ■ (F.59) 


Finally we can read from the previous formula the two point correlator which 
reads 




Ql 


Ql{2)+Ql{l) QiQ 2((2)-(1)) 
Qig2((2)-(1)) Q\{l) + Ql{2) 


(F.60) 


Using the results of section ^ we can hnd the conductivities for each Ai and 
A 2 and then substitute the results in the previous formula to extract the desired 
conductivities associated to Oi and 02 - We obtain 
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(F.61) 

(F.62) 

(F.63) 


It is important to emphasize that we are studying the IR behavior for the con¬ 
ductivities (a; -C 1), therefore, if the constraints on z, 6 allow a negative value for mi 
this term will dominate in the frequency dependence for sufficiently small u. Notice 
also that this computation has been done at zero temperature. On the other hand, at 
hnite temperature we would expect this scaling, to appear as an intermediate scaling 
if T <C a; <C A, with A the interpolation scale for which the IR geometry flow to the 
conformal UV AdS geometry. 

In figure we show the regions of negative mi and m 2 for d = 2,3. We observe 
that for both values of d there is an allowed region of negative mi whereas that 
m 2 has to be always positive, in consequence we conclude that for small enough 
frequencies the conductivities will behave as 


^®The quantities (1) and (2) are only functions of the frequency and the parameters of the system. 
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Figure 14: Contour plots to illustrate the region in the parameter space where the ex¬ 
ponents nil (rainbow coloured) and m 2 (scales of brown) take negative values for the two 
charges model. Left: Conductivity for d = 2. Right: Conductivity for d = 3. The allowed 
values for the parameters are bounded by the gray mesh. The negative values for m 2 are 
outside the permitted region, but there is an allowed region for mi. 
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(F.64) 

(F.65) 

(F.66) 


Notice that mi and m 2 coincide with n and m respectively in the previons snbsection, 
after setting (5 k = 0 in the formnla for n. The reason for this similarity originates 
in the fact that in the previons conhgnration one of the gange helds is charged and 
the charge of the second one is snbleading in the IR, whereas in the present case 
the redehned gange helds correspond to a conhgnration in which the backgronnd has 
charge Qt associated to the hrst gange held, bnt is nentral respect the second one. 
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